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AERONAUTIC  SYMBOLS 
1,  FUNDAMENTAL  AND  DERIVED  UNITS 


Metric 

English 

Symbol 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

Length _ 

1 

meter _  ...  . 

m 

foot  (or  mile) _ 

ft  (or  mi) 

t 

second _ _ 

8 

second  (or  hour) _ 

sec  (or  hr) 

Force _ 

F 

weight  of  1  kilogram - 

kg 

weight  of  1  pound..  . 

lb 

Power  -- 
Speed  —  . 

P 

V 

horsepower  (metric). _ 

fkilometers  per  hour.  . 
(meters  per  second.  - 

kph 

mps 

horsepower.  .  _ 

miles  per  hour.  .  _ 

feet  per  second - - 

Xh 

fps 

2.  GENERAL  SYMBOLS 


Weight = mgr 

Standard  acceleration  of  gravity =9.80065  m/s^ 
or  32.1740  ft/sec^ 

Mass=~ 

.9 

Moment  of  inertia=mt*.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


V  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m“*-s^  at  15°  C 
and  760  mm;  or  0.002378  Ib-ff  ^  sec* 

Specific  weight  of  “standard”  air,  1.2255  kg/m®  or 
0.07651  Ib/cu  ft 


3.  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio,  ^ 

True  air  speed 
Dynamic  pressure, 

f 

Lift,  absolute  coefficient 
Drag,  absolute  coefficient 
Profile  drag,  absolute  coefficient 

D- 

Induced  drag,  absolute  coefficient 
Parasite  drag,  absolute  coefficient  Gi,p=^ 

C 

Cross-wind  force,  absolute  coefficient  Gc=^ 


Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Resultant  moment 
Resultant  angular  velocity 

Reynolds  number,  where  Zis  a  linear  dimen¬ 

sion  (c.g.,  for  an  aiidoil  of  1 .0  ft  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  number  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  number  is  6,865,000) 

Angle  of  attack 

Angle  of  downwash 

Angle  of  attack,  infinite  aspect  ratio 

Angle  of  attack,  induced 

Angle  of  attack,  absolute  (measured  from  zoto- 
lift  position) 

Flight-path  angle 
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SUMMARY 

The  reKiilix  of  a  fheoreficaJ  and  e.nperimenta!  inre.stit/alion  of 
wall  interjerence  for  an  airfoil  i^panninei  a  closed-throat  circular 
wind  tunnel  arc  presented.  Analytical  equations  are  derired 
which  relate  the  characteri.stics  of  aai,  airfoil  in  the  tunnel  at 
subsonic  speeds  with  the  characteri.stics  in.  free  air.  The 
analysis  takes  into  con.sideration  the  effect  of  fluid  compressi¬ 
bility  and  is  based  upon  the  assumption  that  the  chord  of  the 
airfoil  is  .small  as  compared  with  the  diameter  of  the  tunnel. 
The  development  is  restricted  to  an  untwisted,  constant-chord 
airfoil  spanning  the  middle  (f  the  tunnel.  Brief  theoretical 
consideration  is  also  given  to  the  problem  of  choking  at  high 
s])ee(ls.  Results  are  then  presented  of  te.sfs  to  determine  the  low- 
speed,  characteristics  of  an  ArK'Vi.  .fl/liB  airfoil  for  two  chord- 
diameter  ratios.  While,  on  the  basis  if  these  experiments,  no 
appraisal  is  possible  if  the  accuracy  of  the  corrections  at  high 
speeds,  the  data  indicate  that  at  low  Mach  numbers  the  analytical 
reisults  are  valid,  even  for  relatively  large  values  of  the  chord- 
diameter  ratio. 

INTRODUCTION 

The  design  of  nioderu  higli-perfornumee  airi)laiies  requires, 
insofar  as  possibl(',  an  aeeurate  knowledge  of  airfoil  profile 
data  at  Keynolds  and  Maeh  numbers  attained  in  (light. 
Since  the  size  and  ])ower  of  wind  tunnels  are  subj(H't  to  vari¬ 
ous  practical  limitations,  most,  e.xisting  tunnels,  even  if  tlu'y 
can  provide  the  desired  Maeh  numbej',  are  not  capable  of 
attaining  full-scale  Keynolds  numbers  for  all  flight  eonditions. 
To  minimize  this  shortcoming  in  tunnel  tests  of  airfoil 
profiles,  it  is  therefore'  necessary  to  use  moch'ls  having  as 
large  a  chord  us  possible  relative,  to  the  cross-sectional 
dimensions  of  the  tunnel  test  se'ction.  In  order  to  eliminate 
the  efl'ects  of  su]q)ortijig  struts  and  to  ('.xclude  the  ituh'termi- 
nat(‘  tunnel-boundary  interfere-eice  involved  in  the  testing 
of  large-chord  airfoils  of  limited  span,  it  has  bc'come  common 
practice  in  such  tests  to  use  airfoils  which  com])letely  span 
the  test  section.  Eveii  foj’  these*  so-called  “through” 
models,  however,  the'  lunnel-boujidary  intei'ference*  can  still 
be  considerable,  and  accurate  correction  must  be  made  for 
its  efl'ects  if  the  tunnel  data  are  to  be  used  with  confitlence 
in  the  calcnlation  of  free-flight  airplane  characteristics. 

The  tunnel-boundary  interference  for  airfoils  spanning 
wind  tunnels  of  various  types  has  been  the  subject  of  numer¬ 
ous  theoretical  and  exj)erimental  investigations.  The  intei- 
ference  for  rectangulai-  tunnels  having  rigid  walls  normal  to 
the  span  of  the  airfoil  and  either  I'igid  walls  or  free  boundaries 


parallel  to  the  span  has  been  discussed  tlu'oretically  by  sev¬ 
eral  writers.  l^Mr  example,  Lock  (reference  1),  Cllauert  (ref- 
('rence  2),  and  (doldsteii:  (reference  II),  giv<'  the  necessary 
tunnel-wall  eoi-rections  foj-  an  airfoil  spanning  a  rectangular 
tunnel  in  an  incompressible  fluid  ;  while  Cloldstein  and  Young 
(reference  4)  show  how  those  corrections,  as  well  as  those  foi’ 
any  general  case  of  interference  in  an  incompressible  fluid, 
can  be  modified  to  take  account  of  fluid  compressibilit\^ 
Reference  o  gives  tbe  cori'octions  for  tbe  compressible  case 
in  a  closed-throat  rectangular  tunnel,  as  well  as  a  critical 
discussion  of  the  results  of  the  previous  references  and  some 
experimental  data  from  low-speed  tests.  Fage  (reference  G) 
also  presents  experimental  flrag  data  for  several  symmetrical 
bodies  of  various  sizes  in  a  closed-throat  rectangular  tunnel. 
Experimental  and  theoietical  results  for  an  airfoil  spanning 
a  completely  open-thi'oat  rectangular  tunm'l  are  givc'n  by 
8tui)er  (refei'ences  7  aj d  8).  The  ease  of  an  airfoil  spanning 
an  open-throat  circular  tumiel  has  been  the  subject  of  a 
number  of  investigations,  including  theoretical  treatments 
by  Glauert  (referemre  9),  Stiiper  (references  7  and  8),  and 
Squire  (reference  10),  and  expeiimenlal  measurements  by 
Stiiper  (references  7  and  8)  and  Adamson  (reference  11). 
Apparently,  the  case  of  the  closed-throat  circulai'  tunnel  has 
received  no  attention. 

Since  this  last  case  is  often  ('ucountered  in  practice,  an  in¬ 
vestigation  was  matle  of  the  tunnel-wall  interference  at  sub¬ 
sonic  speeds  for  a  wing  spanning  a  closed-throat  (ircular 
tunnel.  The  present  paper  presents  the  results  of  this  in¬ 
vestigation.  In  tile  first  part  of  the  paper,  analytical  equa¬ 
tions  are  tlerived  relating  the  characteristics  of  the  airfoil 
in  the  tunnel  with  those  in  free  air  for  a  compressible  fluid. 
Some  consideration  is  also  given  to  the  phenomenon  of  chok¬ 
ing  which  occurs  at  high  speeds.  In  the  second  part,  the 
validity  of  the  theoretical  results  is  examined  by  the  analysis 
of  experimental  data  for  an  NACA  4412  airfoil  for  two  ratios 
of  airfoil  chord  to  tunnel  diameter.  The  iiivesLigatio  i  is 
restricted  to  Uiitwisted  constant-chord  airfoils  spanning  the 
middle  of  the  tunnel. 

THEORY 

As  in  reference  o,  the  theoretical  development  of  the 
tunnel-wall  corrections  is  divided  conveniently  into  two  gen¬ 
eral  sections.  First,  the  influence  of  the  wall  upon  the  field 
of  flow  at  the  airfoil  in  the  tunnel  is  determined.  Second, 
the  aerodynamic  characteristics  of  the  aii'foil  in  this  field  of 
How  are  related  to  the  corresponding  quantities  in  free  air. 
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Id  this  maniu'i-,  simple  formulas  are  finally  obtained  which 
enable  the  prediction  of  the  free-air  characteristics  when  the 
characteristics  in  the  tunnel  are  known. 

Again,  as  in  reference  5,  the  analysis  is  based  upon  the 
method  of  superposition.  To  this  end,  it  is  assumed  that  the 
airfoil  is  of  small  thickness  and  camber,  so  that  the  induced 
velocity  is  everywhere  small  as  compared  with  th(‘  velocity 
of  the  undisturbed  stream.  With  this  assumption,  the  total 
induced  vidocity  at  any  point  is  the  simple  vectoi'  sum  of  the 
sejiarate  velocities  induced  at  that  point  by  the  interference' 
hetwe'en  the  tunnel  wall  and  tlu'  airfoil  camber,  thickness, 
and  wake.  Thus  the  efl'ects  of  cambi'r,  thickness,  and  wake 
may  ('ach  be  analyze'd  sepai'ately  and  superposed  to  obtain 
th('  desired  I’esult  for  the  complete  airfoil.  As  pointed  out 
in  ref('rence  5,  this  procedure  is  pe'rmissible  even  in  the  com¬ 
pressible  thud  if  the  airfoil  is  of  small  thickness  ami  camber 
as  assumed. 

Ih'fore  proceeding  to  the  actual  development  of  tlu'  iheoi-y, 
it  is  usc'fid  to  contrast  the  pres('nt  problem  with  tin'  probh'ins 
of  through  ail-foils  in  the  various  types  of  rectangular  tunnels 
and  in  the  open-throat  circular  tunnel.  In  the  case  of  an 
airfoil  spanning  a  rectangular  tunnel  having  rigid  walls 
normal  to  the  span  of  the  airfoil,  tlm  problem  is  ri'latively 
simph'.  If  the  effect  of  the  boundary  layer  along  tin'  lunnel 
walls  is  neglected,  the  flow  is  sensibly  the  same  in  all  plain's 
normal  to  the  span;  that  is,  there  is  clearly  no  spanwise 
variation  in  lift.  The  air  flow  is  thus  essentially  two-dimen¬ 
sional,  and  the  interference  problems  of  camber,  thickness, 
and  wake  can  be  analyzed  by  tin'  customary  means  of  a 
systc'in  of  images  with  axes  [laralli'l  to  the  span  of  tin'  airfoil 
(references  1,  2,  3,  and  5).  This  is  true  whetln'r  tin'  tunnel 
boundaries  parallel  to  the  span  ari'  fixc'd  or  fri'i'.  In  this 
manner,  tunnel-boundary  corrections  can  be  derived  for 
airfoils  of  moderately  large  chord  as  compared  with  the 
height  of  the  tunnel  test  section. 

In  tin'  case  of  an  airfoil  spanning  a  coin])let('ly  fri'i'  ji't, 
wlu'tln'r  rectangular  or  circular  in  sc'ction,  the  lift  in'c<'ssarily 
falls  to  zero  at  the  boundary  of  tin'  jet.  1  hci'e  thus  ('xists 
in  this  case  a  pronounced  spanwise  vai'iation  in  lift  and  an 
attc'udant  syst('m  of  trailing  vortices.  In  the  existing  treat- 
nu'nts  of  the  probh'm,  oidy  the  intcj'fi'rence  between  these 
trailing  vorticc'S  and  the  jet  boundaric'S  is  C'onsidered,  the 
interference  effects  associat('d  with  the  chordwisc  distribu¬ 
tion  of  bound  vortices  and  with  tlu'  airfoil  thickness  and 
wak('  being  completely  negh'ctcd.  This  procedure  implies 
tin'  assumption  that  the  ('hord  of  tin'  airfoil  is  V('ry  small 
i-('lativ('  to  tin'  dimensions  ol  the  jet.  In  this  manin'r,  tin' 
problem  is  reduced  to  a  limiting  case  of  tin'  usual  jjroblem 
of  an  airfoil  pai  tially  spanning  the  jet,  and,  as  in  this  latter 
case,  the  component  of  downwash  induced  at  the  airfoil  by 
the  interference  between  the  walla  and  the  trailing  vortices 
is  one-half  as  great  as  the  corresponding  component  an 
infinitt'  distance  downstr<'am.  The  theoretical  ih'tcj-mina- 
tion  of  the  wall  interference  may  thus  be  treated  as  a  ju’oh- 
h'ni  of  two-dimensional  flow  in  a  plane  normal  to  the  axis  ol 
the  tunnel  infinitely  fai-  behind  the  airfoil.  The  boundai-y 
cimditions  for  ('itlu'r  tlu'  rc'ctangulai-  or  circular  j('t  ai’c  then 
readily  satisfied  by  tin'  introduction  of  a  suitable  system  of 
image  vortices  with  axes  parallel  to  the  axis  of  tin'  tunm'l 


(references  7,  8,  9,  and  10).  This  method  of  analysis,  how¬ 
ever,  is  inadequate  if  the  chord  of  the  airfoil  is  even  moder¬ 
ately  lai'ge  as  compai-ed  with  the  dimensions  of  the  jet. 

The  case  of  the  airfoil  spanning  a  closed-throat  circular 
tunnel  is  more  complex  than  eitlu'r  of  the  foregoing  problems. 
Unlike  the  condition  prevailing  in  the  free  jet,  the  lift  in 
this  case  need  not  fall  to  zero  at  the  boundary — that  is,  at 
the  tunnel  wall — so  that  the  spanwise  variation  in  lift  is 
not  necessarily  large.  In  fact,  as  will  be  seen,  the  lift  is 
constant  across  the  span  of  the  airfoil,  and  no  sysb'in  of 
trailing  vortices  exists.  The  assumption  of  a  very  small 
chord  and  the  consequent  reduction  of  the  problem  to  a  case 
of  two-dimensional  flow  in  a  plane  infinitely  far  dow  iistri'am 
is  thus  wdthout  meaning.  On  the  other  hand, an  analysis 
for  airfoils  of  moderately  large  chord  in  the  mannei-  em- 
ployetl  in  the  case  of  the  rectangular  tunnel  with  rigid  side 
walls  is  not  possible.  In  the  closed-throat  circular  tunni'l 
th('  flow-  in  all  planes  normal  to  the  span  of  the  airfoil  is  not 
the  same,  so  that  the  eftect  of  the  bound  vortices,  and  of  tlu' 
airfoil  thickness  and  wake  as  well,  cannot  be  treated  as  a 
probh'in  in  t\sn-dimensional  flowu  Furthermore,  the  bound¬ 
ary  conditions  at  the  tunnel  wall  cannot  be  satisfied  for  the 
actual  three-dimensional  problem  by  any  knowm  system 
of  images.  The  solution  of  the  problem  for  the  closed-throat 
circular  tunnel  thus  requires  an  analysis  entirely  different 
from  those  employed  in  the  previous  instances. 

INFLUENCE  OF  TUNNEL  WALL  UPON  FIELD  OF  FLOW  AT  AIRFOIL 

An  approach  to  the  problem  of  the  aii-foil  spanning  a 
('losed-throat  circular  tunnel  is  afl’oi-ded  by  the  work  of  von 
Karman  and  Burgers  in  ri'ferenci'  12  (pj).  266  to  273),  wdiei-c' 
the  velocity  potential  at  an  arbitrary  point  in  the  tunnel  is 
determined  for  a  U-shape  vortex  of  infinitesimal  span  in  an 
incompi essible  fluid. 
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A  syslcin  of  roct angular  coordinates  x,  y,  z  is  introduced 
as  shown  in  figuri'  1 .  The  x-axis  is  taken,  on  the  center  liiu' 
of  the  tunnel  with  its  ]iositive  direction  downstream.  Tlu' 
2-axis  is  positive  downward,  and  the  )/-axis  jiositivo  to  the 
h'ft  for  an  observer  looking  against  the  direction  of  flow. 
An  alternativi'  system  of  cylindrical  coordinates  x,  co,  6  is 
tlelined  by  tlu'  relations 


y=<jo  cos  d 
z=w  sin  6 


(I) 


The  positive'  diri'ctiou  of  circulation  is  fleliiu'd  so  that  a 
voi'tex  with  positive  circulation  exerts  a  force  on  the  fluid  in 
the  direction  of  the  ]jositiv('  2-axis.  In  other  words,  the  lift 
force  experienced  by  a  positive  voi'tex  is  in  the  negative  z 
direction.  The  velocity  of  the  fluid  in  the  undisturbed 
stri'am  is  denotc'd  by  Y'  and  the  radius  of  the  tunnel  by  r. 
Other  symbols  are  defined  as  introduced  in  tbe  te.xt.  A  list 
of  the  more  important  symbols  and  their  definitions  is  given 
in  Appendix  C. 

(Auisider  now  a  U-sliape  vorti'x  of  iufinitisimal  span  dr\ 
liarallel  to  tlu'  )/-axis  and  situati'd  in  the  ?/2-j)lane  at  the  point 
n  —  cos  00,  f  =  sin  0().  If  th('  strength  of  the  vortc'x  is 
denoted  by  F'  the  velocity  potc'iitial  in  the  closed  tunnel  at 
the  points  x,  w,  6  is  given  by  von  Karman  and  Burgei's,  for 
negative  vahu's  of  x,  as 


i-rfW-oa 

TT  J(l  df 
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(It  should  be  noted  that  the  quantity  ^  appearing  in  these 
equations  is  merely  a  variable  of  integration  and  has  no  physi¬ 
cal  significance.)  The  quantity  is  a  Bessel  function 

of  the  first  kind  of  the  order  m.  The  summation  with  respect 
to  rn  extends  over  all  the  positive  integers  and  includes 
m  =  the  prime  added  to  the  summation  sign  indicates  that 
a  factor  h  must  be  inserted  before  the  term  corresponding  to 
w=0.  The  summation  with  respect  to  .s'  for  every  m 
extends  over  all  positive;  roots  of  the  equation 

t/m'(k.iC)  =  0  (4) 

where  is  the  derivative  of  the  function  with 

respect  to  its  argument.  The  notation  used  throughout  this 
paper  for  the  Bessel  functions  is  that  of  Watson  (reference 
Id),  which  is  the  same  as  that  of  the  Smithsonian  table's 
(reference  14). 

By  differentiating  Q  with  respect  to  f  and  then  intiigrating 
with  respect  to  as  indicated  in  eejuation  (2),  the  velocity 
potential  becomes  finally 
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As  pointed  out,  this  expression  ajiplies  only  at  negative  valiu's 
of  X.  As  will  be  seen  later,  the  necessary  results  for  positive 
values  of  x  can  be  derived  from  considerations  of  symmetry. 

By  means  of  (filiation  (5),  it  is  possible  to  evaluati'  the 
wall  interference  associated  with  both  airfoil  camber  and 
thickness  for  the  case  of  the  incompressible  fluid.  Tlu'sc 
results  can  then  be  modified  for  the  effect  of  fluid  com¬ 
pressibility  by  the  methods  of  reference  4.  It  is  found 
finally  that,  for  a  closed-throat  circular  tunnel,  the  effects  of 
interference  between  the  walls  and  the  airfoil  camber  are 
identical  with  the  corresponding  interference  effects  for  the 
same  airfoil  spanning  a  closed-throat  rectangular  tunnel,  tlie 
height  of  which  bears  a  known  relation  to  the  diameter  of 
the  circular  tunni'l.  A  similar  conclusion  is  obtained  r(;gard- 
ing  tbe  ('fl'ects  of  interference  between  tlu'  walls  and  tlu'  air¬ 
foil  thickness,  except  for  a  numerical  dift'erence  in  the 
relation  between  the  diamctei'  of  the  given  circular  tunnel 
and  the  height  of  the  equivalent  rectangular  tunnel.  The 
interference  effects  associated  with  the  wake  of  the  airfoil 
are  not  analyzed  in  detail,  but  their  magnitude  can  be 
estimated  with  reasonable  accuracy  by  comparison  with  the 
results  for  the  tliickness  effect.  In  order  to  simplify  the 
complex  mathematics  of  the  problem,  the  interference  (flfects 
are  calculated  only  for  the,  section  of  the  airfoil  at  the  center 
line  of  the  tunnel.  As  will  bo  seen  later,  however,  experi¬ 
mental  data  indicate  that  the  results  are  applicable  at  any 
spanwise  station. 

Camber  effect. — To  analyze  the  effect  of  the  interference 
between  the  tunnel  walls  and  the  airfoil  camber,  the  thi(;kness 
and  wake  of  the  airfoil  are  considered  to  be  removed  and  the 
airfoil  reduco'd  to  its  mean  camber  line.  The  resulting  infini¬ 
tesimally  thin  airfoil  may  then  be  replaced  by  a  sheet  of 
continuously  distributed,  bound  vortices  which,  in  the  general 
three-dimensional  casi',  consist  of  both  spanwise  and  chord- 
wise  vortices.  The  velocity  induced  at  any  given  point  on 
the  camber  line  is  then  obtained  by  integration  over  the  entire 
vortex  sheet.  As  in  all  thin-airfoil  theory  the  distribution 
of  bound  vorticity  must  be  such  that  the  resultant  of  this 
induced  velocity  aiul  the  fi'ee-stream  velocity  is  tangential 
to  the  camber  line  at  all  points.  As  will  be  seen,  however, 
the  actual  theoretical  determination  of  the  distribution  of 
vorticity  is  not  necessary  in  this  case. 

In  calculating  the  velocity  field  of  the  vortex  system,  it  is 
assumed  that  the  bound  vorticity  is  distributed  in  the  middh; 
plane  of  the  tunnel — that  is,  in  the  .ry-plane  — rather  than 
along  tlie  camber  line  and  that  the  induced  velocity  at  any 
point  on  the  camber  line  is  the  same  as  the  induced  velocity 
at  the  corresponding  point  in  the  xy-plane.  From  equation 
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(5),  tlip  Toloritv  poti'iitial  ai  any  point  .r,  co,  d  for  a  vortex 
eh'ment  on  tlie  y-nx\s  at  tin'  ])oint  v  (8n=i),  ^0=17)  is 
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The  term  for  m  =  0  disappears  b\'  virtue  of  the  faetor  m  in 
the  nnmei-ator  of  the  general  term.  The  vertieal  indiieed 
vi'loeity  r/  in  tin'  incompressible  Hnid  is  then 
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The  complicated  double  series  in  this  equation  can  be 
reduced  to  a  single  seric's  and  the  mathematics  of  the  problem 
greatly  simplified  by  limiting  the  discussion  to  the  chord- 
wise.  section  of  the  airfoil  at  the  center  line  of  the  tunnel 
(?/==0.  2  =  0),  From  the  known  relations  for  Bessel  functions 
(cf.  reference  1 3),  for  7/  =  0 


J,(V,S)_|,2 
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Thus,  at  points  on  the  .r-nxis, 
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where  the  summation  with  respect  to  .9  extemis  over  all  the 
positive  roots  of  the  equation 

-/i'(Xpd  =  0  (9) 

From  Bessel’s  differential  equation 

J/'(Kr)  --(l-^l,,2)-/d>»'-}  (lO) 


where  the  double  prime  denotes  the  second  derivative  of  the 
Bessel  function  with  respect  to  its  argument.  Equation  (8) 
can  thus  be  written 

/ 1  _  dv  c  ^  tJi  {\sv)  tilt 

"  27rr77  Xsrt/i(X,,r)  Ji''(X,,r)  ■ 


As  mentioned,  this  equation  is  valid  only  for  negative  values 
of  X. 

It  is  apparent  that  the  series  in  equation  (11)  is  rajiidly 
convergent  foi-  large  negative  values  of  x,  but  that  the  con- 
vei’gence,  is  slow  for  small  negative  values.  Since  in  iJie 
evaluation  of  the  velocity  induced  by  the  voi’tex  sheet  it 


is  the  small  values  of  x  which  are  of  primaiy  importance, 
equation  (11)  cannot  be  applied  directly  in  the,  present  case. 
It  is  possible,  however,  b\'  means  of  a  method  demonstrated 
by  Watson  (reference  15),  to  exjirc'ss  the  series  of  this  equa¬ 
tion  as  a  combination  of  ('lementary  functions  and  a  series 
of  ascending  powers  of  x  and  rj.  The  resulting  sc'ries  is 
readily  applicable  to  the  jnesent  problem. 

The  detaileil  procedure  for  the  transformation  of  the  series 
of  e(|uation  (11)  is  given  in  Appendix  A.  By  application 
of  the  final  I'csult,  ('([nation  (II)  may  be  written. 
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'I’he  double  summation  extends  over  all  inti'gral  values  of  k 
and  p  from  zero  to  [lositive  infinity.  ’Phe  numerical  coeffi¬ 
cient  M%((+/<+i)  =  /c^2/  is  given  by  the  integral 
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Here  I,{t)  is  a  modified  Bessel  function  of  the  first  kind  of 
order  unity,  and  i/(0  denotes  the  derivative  of  /i(f)  with 
respect  to  its  argum('nt.  The  numerical  values  of  ,u'2/  for 
/=1,  2,  3,  4,  arc  evaluated  by  means  of  a  series  expansion 
in  Appendix  A. 

It  is  readily  shown  that  the  first  term  on  the  right-hand 
side  of  equation  (12)  agrees  with  the  induced  velocity  com¬ 
puted  for  x=0  by  the  more  elementary  theory  of  tunnel-wall 
interference  which  considers  only  the  eft'ccts  of  the  trailing 
vorli(,‘es  and  their  images.  To  this  end,  consider  the  two- 
dimensional  flow  in  a  plane  normal  to  the  axis  of  the  tunnel 
an  infinite,  distance  downstream  (fig.  2).  'Fhe  theory  states 
t  hat  the  induced  v(‘locity  at  a  given  point  (?/,  2)  in  this  plane  is 
twice  as  great  as  the  ind  i  iced  velocity  at  th  e  corresponding  point 
in  the  plane  r=0  (cf.  reference  12,  p.  260).  In  the  plane 
X--  00,  the  trailing  vortices  of  the  U-shape  vortex  previously 
considered  constitute,  a  vortex  pair  having  an  infinitesimal 
spacing  drj  and  situated  at  the  point  y  =  ri,  2  =  0.  The  circu¬ 
lation  of  each  voi-tex  of  the  pair  is  F'  and  is  directed  as  shown 
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KiBf'RK  2.  Section  through  luiiiiol  at  infinity  downsfream. 
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in  figure  2.  The  boundary  condition  that  tliero  sliall  be  no 
flow  normal  lo  the  wall  of  the  tunnel  can  be  satisfied  by  the 
introdiietion  at  llu'  ])oint  z=0,  ?/=/'-/■()  of  an  image  vortex 
])air  with  a  spacing  r’d-qlr]'^  and  with  the  eirculatioti  of  the 
vortices  directed  as  indicated.  The  vertical  velocity  induced 
at  the  tunnel  cento'r  by  the  trailing  voi'tex  pair  is 


and  the  vertical  velocity  induced  at  llu'  same  ])oint  by  llic 
image  vortex  pair  is 


T'dr, 


'I'lie  total  vei'tical  velocity  at  the  center  of  the  tunnel  at 
x=  <x>  is  then  the  sum  of  these  two  velocities;  that  is, 


T'dri(r~d~‘n^) 
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The  vertical  velocity  at  the  center  of  the  tunnel  at  x^  i)  is 
one-half  of  this  value,  or 
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This  value  agrees  with  the  result  of  equation  (12)  for  the 
special  case  x=0.  Thus,  the  first  term  on  the  right-hand 
side  of  equation  (12)  represents  the  vertical  induced  velocity 
on  the  center  line  of  the  tunnel  at  x=0  and  is  attributable 
entirely  to  the  trailing  vortices  and  to  the  interference 
between  these  vortices  and  the  tunnel  w^alls.  The  remaining 
terms  represent  the  variation  in  induced  velocity  duo  to  a 
displacement  a  distance  x  upstream  from  the  origin.  These 
terms  arise  both  from  a  change  in  the  effect  of  the  trailing 
vortices  and  their  wall  interference  and  from  the  now-active 
effect  of  the  transverse  bound  vortex  and  its  interference 
with  the  tunnel  walls. 

Although  equation  (12)  was  deduced  for  negative  values 
of  X,  it  can  be  showui  that  it  is  applicable  to  positive  values 
of  X  as  well.  According  to  von  Karman  and  Burgeis  (ref¬ 
erence  12,  p.  267),  the  vortical  induced  velocity  at  — r  is 
relat('d  to  the  corresponding  velocity  at  -hr  by  the  equation 

r/(— r)  =  r/(oc)— c/(  +  r) 

By  virtue  of  this  relation,  together  wdth  the  fact  that 
r/(0)  p/(  00  ),  it  follow's  that 


r/ ( +  x)  -  r/  (0)  =  -  [r/  ( -  x)  -  r/  (0 )  ]  (16) 


That  is,  the  difference  between  the  induced  velocity  at  a 
given  station  x  and  the  induced  velocity  at  r=0  must  be. 
an  odd  function  of  x.  The  terms  containing  x  in  ecpiation 
(12),  which  were;  derivt'd  to  represent  this  difference  for 
negative  values  of  the  variable,  are  seen  to  constitute  pre¬ 
cisely  such  a  function.  Thus  the  expansion  of  equation  (12) 
is  valitl  for  positive  as  well  as  m'gative  values  of  x. 


The  vortex  sheet  which  repres('nts  the  entire  airfoil  can 
now"  be  built  up  by  the  superposition  of  U -shape  vortices  in 
the  x?/-plane,  and  the  total  induced  velocity  found  by  inte¬ 
gration  of  ecpiation  (12)  oviu'  the  entire  system.  The  leading 
('dge  of  the  airfoil  is  placed  on  the  ?/-axis  as  shown  in  figure  3; 
the  ti-ailing  edge  then  lies  at  x=c,  where  c  is  the  chord  of  the 
airfoil.  The  circulation  of  an  ehmentary  vortex  having  an 
infinitt'simal  span  r/17  and  situated  at  the  point  x=?,  y=n  if’ 
taken  to  b(>  (dr'ld^)d^,  wIktc  (c/r'/d?)  is  the  vorticity  per 
unit  length  of  tlii'  chordwise  section  at  the  station  y  =  rj. 
Th(>  vertical  velocity  induced  at  the  chordwdse  station  x  on 
th('  centc'r  line  of  the  tunnel  by  a  single  elementary  vorte.x  is 
given  by  equation  (12)  if  x  and  T'  are  replaced  by  (x— ^) 
and  {dr'ld^)di,  respectively.  The  total  vertical  vdocity 
induced  by  the  complete  airfoil  is  then  given  by  the  double 
integral 
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The  integration  of  ecpiation  (17)  requires  a  knowledge  of 
(dr'/ds)  as  a  function  of  rj  and  Theoretically,  {dV'jd^) 
coidd  be  determined  from  the  requirement  that  the  induced 
vertical  velocity  at  every  point  on  the  camber  line  must 
be  such  that  the  resultant  of  this  velocity  and  the  lrc(‘- 
stream  velocity  is  tangential  to  the  camber  line.  This 
method  of  procedure  leads,  however,  to  a  complicated  double 
integral  equation,  the  solution  of  wliich  does  not  appear 
feasible.  Some  assumption  concerning  the  distribution  of 
vorticity  must  therefore  be  made  if  the  problem  is  to  be 
solvixl . 
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To  aid  ill  llic  clioice  of  a  suitable  assumption,  experiments 
were  carried  out  to  determine  tin'  jiressure  distribution,  both 
ehoi'dwise  and  spanwise,  over  an  airfoil  spanning  a  closed- 
throat  circular  tunnel.  The  airfoil  used  in  the  experiments, 
which  ar(>  descrihed  in  detail  later  in  the  report,  had  an 
NACA  4412  section  and  was  untwisted  and  of  constant  chord. 
'Phe  results  of  these  experiemnts  reveal  that  for  such  an 
ai'rangement  the  lift  is  sensibly  uniform  across  the  span  for 
angles  of  attack  lielow  the  stall.  This  fact  is  illustrated  in 
figures  (>  and  7,  which  show  the  experimental  spanwise 
lift  distribution  for  the  airfoil  at  various  angles  of  attack  in 
wind  tunnels  afToi’ding  chord-diameter  ratios  of  0.3,57  ami 
0.625.  These  results  were  at  first  regarded  as  rather  sur¬ 
prising.  Later,  however,  it  was  realized  that  they  are  only 
what  might  logically  he  expected  from  general  considei-ations 
of  tlie  conditions  of  flow  in  a  closed-throat  tunnel.  A 
demonstration  of  this  fact  is  given  in  Appendix  B,  in  which 
it  is  showji  that  tlu'  lift  distribution  is  uniform  across  a.n 
untwisted,  constant-chord  airfoil  spanning  any  clos('d-f hroat 
wind  tunnel,  irrespective  of  the  cross-sectional  shape  of  the 
tunnel.  Detailed  examination  of  the  ])r(>ssure  distributions 
from  which  the  results  of  figures  6  and  7  were  obtained 
reveal  further  that  at  a  givem  angle  of  attack  tin'  chordwise. 
pressure  distribution  is  sensibly  the  same  for  all  spanwise 
stations  on  the  airfoil;  that  is,  the  lift  per  unit  chord  at  any 
given  chordwise  station  is  constant  across  the  span.  It  is 
to  be  expeched  that  this  result,  though  obtained  for  a  par¬ 
ticular  airfoil,  will  he  equally  true  for  any  ordinary  eamber- 
line  shape.  Thus  it  is  reasonable  to  assume  that  the  dis¬ 


tribution  of  bound  vorticity  is  not  a  function  of  the  spanwise 
position  on  the  airfoil;  that  is.  Wr'/c/|)  is  independent  of  rj. 

On  the  basis  of  this  assumption,  e<(uabion  (17)  may  be 
written 
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and  th(>  integration  carried  out  with  respect  to  r;.  The  first 
two  terms  of  the  integi'and,  however,  become  infinite  at  th(' 
])oint  jj=0.  These  singularities,  which  are  due  to  the  effects 
of  the  vortices  trailing  from  the  vmrt('x  elements  on  the' 
T-axis,  i-equire  that  special  carc'  be  taken  in  the  integration. 
Tlu'  evaluation  of  the  integral  must  bc'  carried  out  from  —  r 
to  — 6  and  from  e  to  +r,  and  to  the  resulting  function 
must  b(‘  added  the  effects  of  the  trailing  vortex  pairs  of  span 
2  6  which  straddle  the  /-axis.  The  limit  of  this  sum  must 
then  l)e  taken  as  e  terids  to  zei’o.  The  vertical  velocity 
induced  at  the  point  /  on  the  /-axis  by  tin'  vortices  trailing 
from  a  vortex  element  of  span  2e  .symmetrically  placed  at 
,'/=0  is 
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Since  the  first  two  terms  of  equation  (18)  contain  only 
s('con<l-order  powers  of  rj,  th('  integrals  from  —  r  to  —  e  and 
from  -1-6  to  +r  will  he  equal.  The  integral  of  these  tw’o 
terms  with  respect  to  17  thus  becomes  finally 
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The  integration  with  inspect  to  r]  of  the  double  series  in 
equation  (18)  presents  no  difficulty.  The  expression  for 
r/  thus  becomes  after  iiitegration 


h'lit  to  assuming  that  powers  of  the  chord-diameter  ratio 
(c/d)  higher  than  tlie  second  may  be  neglected  in  the  final 
('(piations  for  the  tunnel-wall  corrections.  The  approxi¬ 
mation  is  accom])]ished  by  expanding  the  first  term  of  the 
integraiul  in  ascending  powers  of  (/— ^)/r  and  discarding  all 
terms  containing  powers  higher  than  the  first  and  by  retain¬ 
ing  only  the  p  =  Q  terms  of  the  double  series.  This  gives  for 
the  induced  velocity 
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To]-  constant  spanwdse  circulation,  the  ti'ailing  vortices 
finally  disappear  in  the  integration  with  respect  to  77.  The 
integrand  of  equation  (19)  thus  represents  the  increment  of 
vertical  velocity  induced  by  an  elementary  vortex  of  constant 
circulation  completely  spanning  the  tunnel. 

It  will  now  b('  assumed  that  the  chord  of  the  airfoil  is 
small  enough  as  compared  with  the  dimensions  of  the  wind 
tunnel  that  powers  of  (/— |)/r  greater  than  the  first  may  be 
neglected  in  tlu'  integrand  of  equation  (19).  This  is  ('(piiva- 
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JVv  siibstituliii^^  tlio  numerical  value's  foi'  tlu'  coelficienls 
m'2(/,j.i)  from  ('(|uatioTis  (A20)  of  Ajjjx'udix  A,  this  equation 
Jiiay  he  written  to  an  aceui’aey  of  tlirc'c  sig7nficaiit  figures  as 


0.579 


(it 


(20) 


'riie  foi’egoing  result,  which  was  derived  by  assuming  the 
fluid  to  b('  incompressible,  can  be  modified  for  the  effect  of 
comjn-essibility  by  the  methods  of  Goldsb'in  and  Young. 
'Phe  modification  is  most  readily  performed  hy  means  of  the 
so-called  “Method  11’’  (refei'once  4,  pp.  5-0),  which  compai'es 
the  compressihle  and  incompressible  flows  for  equal  values  of 
circulation.  If  the  Mach  numbi'r  of  the  compressible  flow 
at  tlu'  position  of  the  airfoil  is  dt'nofed  hy  d/,  it  is  readily 
shown  on  the  basis  of  this  method  that  foi'  a  given  distribu¬ 
tion  of  vorticity  tin*  vc'rtical  velocity  induced  in  a  cojupres- 
siblc  fluid  at  any  point  on  the  center  liiu'  of  a  tunnel  of  radius 
r  is  yl— d/-  tinn's  the  corresponding  velocity  at  tlu'  sain*' 
])()int  in  an  incompi'cssible  fluid  in  a  tunnel  of  radius 
/'y  l  — d/-.  Thus,  from  equation  (20),  the  vertical  velocity 
i\'  in  a  coinpressibh'  fluid  in  the  actual  tunnel  of  radius  r  is 


\  l-d/2 

2ir 


.r(' 


dV 


hr  1 


di  (2  1 ) 


Th('  first  term  of  this  equation  I'opresents  the  vertical  velocity 
that  would  b('  induced  by  a  vorte.x  she('t  of  infinite  sjian  in  an 
unlimited  fluid  field.  The  second  term  thus  represents  the 
interference  effect  of  the  tunnel  wall. 

Equation  (21)  may  be  compared  with  the  corresponding 
result  from  reference  5,  which  discusses  the  wall  inti'rference 
for  an  airfoil  in  a  closed-tbroat  two-dimensional-flow  wind 
tunnel.  After  alteration  to  conform  with  (he  notation  and 
sign  conventions  of  the  present  paper,  equation  (41)  of  refor- 
('iice  5  gives  for  the  vertical  velocity  at  the  camber  line  of  an 
inlinitesimally  thin  airfoil  mounted  on  the  ci'iiter  line  of  a 
two-dimensional-flow  tunnel  of  height  li 


r  /dv 


-TT  Jn 


1 


d^J[_x~i  0/h(l 


Comparison  of  eciuations  (21)  and  (22)  shows  that  an  intini- 
tesimally  thin  airfoil  spanning  a  closed-throat  circular  tunnel 
of  radius  r  expei'iences  at  its  midspan  section  the  same  inter- 
f('r('nce  as  would  be  experienced  by  the  same  airfoil  in  a  closed- 
tlu’oat  two-dimensional-flow  tunnel  of  height 


//,= 


T 

t)(0.579) 


r=^  1 .080/' 


or,  in  terms  of  the  tuiuu'l  diameter  d, 


Ai  =  0.84:fd 


(T!) 


'Phis  result  makes  the  later  determination  of  the  interference 
cori’ections  for  the  circular  tunnel  very  simple,  since  the  cor- 
jx'ctions  for  the  rectangular  tunnel  are  already  known. 

It  is  readily  shoevn  by  means  of  equation  (fi)  (hat  the  vortc'x 
system  which  represents  the  infinitesimally  thin  airfoil  in¬ 
duces  no  axial  velocity  at  any  point  in  the  .r?/-plane.  It 
follows  that  airfoil  camber  has  no  effect  upon  (he  axial  ve¬ 
locity  or  pressure  gradient  at  the,  i^osition  of  (he  model. 
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Thickness  effect. — 'Phe  interference  effc'cts  associated  with 
airfoil  thickness  can  b('  found  hy  reducing  the  given  airfoil 
to  its  base  profile  and  analyzing  tlu'  interfer('uc('  b('twe('n  the 
UinrK'l  wall  and  this  profile.  'Phe  base  profile'  is  defnu'd  as 
the  profile  the  airfoil  would  have'  if  tin;  cambc'i-  were'  r('me)ve'd 
and  the  resulting  airfoil  place'el  at  zero  angle  of  attack.  If  it 
is  assumed  that  no  wake  is  pi'e'sent,  the  inteilei-ence'  hctwe'e'ii 
the  tiiiuiel  wall  and  this  symmetrical  airfoil  can  be  founel  by 
applying  the  results  of  ('(luation  (5)  to  Lock’s  me'thoel  of 
analysis  of  the  interfere'nce'  on  a  symmetrical  boelv  in  twee- 
dimejisional  incompressible'  flow.  (Lock’s  oi'iginal  analysis 
appears  iji  referencee  1 ;  an  alternative  explanation  of  the' 
method  is  given  by  Glaiu'rt  in  refeu’cnce  2,  pi).  52-57.) 

Lock’s  method  of  analysis,  which  assunu's  that  the  chore! 
of  the  airfoil  is  small  as  compared  with  the'  elirnensions  of  the 
tunnel,  consists  e'ssentially  in  replacing  the  given  symme'trie'al 
airfoil  by  an  eepiivalent  two-dimcuisional  source-sink  eloubh't 
anel  calculating  thee  interference  betwee'Ji  this  eloublet  aiiel 
the  tunne'l  boiindarie's.  The  stremgth  of  the  doublet  in  any 
given  case  is  proportioneel  so  that  it  ineluces  at  a  considerable 
elistance  frenn  itse'lf  in  free'  air  a  velocity  e'qual  to  the  veleecity 
ineluced  at  the  same'  point  by  the  or-iginal  airfoil.  In  the' 
two-elimensiona)  case,  the  interference  flow  at  the  position 
of  the  airfoil  is  them  r('ae;lily  found  by  introducing  an  infinite' 
se'ric'S  of  image's  of  the'  eloubh't  sue'h  as  to  satisfy  the'  conelitie)n 
that  tlu're  shall  he'  no  flow  normal  to  the  tmuie'l  leounehirie's 
anel  calculating  the  ve'loe'ity  indueieed  at  the  airfoil  by  this 
system  of  images.  For  aei  airfoil  spanning  a  closed-throat 
I'ectangular  tunne'l  at  miel-height,  the  Jiet  result  of  the  wall 
interference  for  the  ine;om press! ble  case  is  to  increase  the^ 
e'ffe'ctive  axial  veloe.'ity  at  the  position  of  the!  airloil  by  the 
amount 


wheix'  ju  is  (he  strength  eef  the'  doulflet  use'd  (e)  repre'se'iit  the' 
;iirfe)il.  It  is  shown  in  re'fei'e'iices  4  iind  5  that  thei  e'ffee't  of 
fluid  compre'ssibility  is  to  ine'i'ease  (his  inte'rfere'nce'  ve'locity 
by  the  factor  1 /|1  —  (d/')^]^'^,  where  Af'  is  the'  Mach  nunda'i' 
e)f  (he  unelisturbe'd  stream  in  the  tunne'l.  Thus,  in  the'  coni- 
])ressible  e'ase, 

A  V'';'--  . -  JhT-.  ..  -  pein 

'PIh'  problem  eef  the  symmetrical  airfoil  in  a  cleese'el-threeat 
e'ircular  (iimiel  e-au  also  be  solved  by  re'plae'ing  the'  lurfeeil  by 
iui  eepnvalent  doublet  spanning  the  tunnel.  In  this  case', 
though,  the  interfei'ence'  for  the  doublet  cannot  be  founel  by 
(he  methoel  of  images.  If  the  eloublet  used  is  compose'd, 
however,  of  two  vortices  in  a  plane  normal  to  the  steam  in- 
ste!ad  of  the  customary  source  and  sink  in  line  with  tbe 
stre'am,  the  interference  vedocity  can  be  calculated  by  nu'ans 
of  ccpiation  (5).  Siiu'c  the  velocity  fie'lds  of  the  two  ty|)e's 
of  doublets  arc  identical,  the  interference'  calculateel  by  nu'ans 
of  the  xmrtex  doubh't  is  the  same  as  that  which  woedel  be' 
obtained  if  the  source',-sink  eloublet  were  useel. 

Consieh'r  a  vortex  e'lement  of  circulation  P'  anel  span  dr] 
at  the  point  a)„,  do  in  the  ?/2-plane  (fig.  1).  Ej'om  equation  (d), 
the  streamwise  ve'loe'ity  o/  indiiceel  at  any  peeint  x,  w,  0  iij)- 
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Stream  from  the  origin  hy  this  element  and  its  accompanying 
trailing  vortices  is 


C)(t> 

d.r 


t  f/rj 

’r/'"  »i=(i 


r^*V,„(X.co) 


X 


/I  ■  //I  /I  \  I 

?fi  cos  sin  /fi(0  —  0(i)  -  -  + 

“o 


Xs  sin  (t,|COs  ni(6*— (/„)  ./, 


,'(X,,w,i) 


(25) 


At  a  jjoint  in  tin'  middle  |)lanc  of  the  tunm'l  (:ny-[)lanc), 
0  =  0,  co==7/,  and  tlu'  velocity  is 


TT/ 


c^;L/,„(X,y)  ^ 

(l-jJ”^)x,J„d(X»r) 


[■ 


>fl  cos  0,1  sin  me,,  ^  j.ijj  ,.„g  JJ  (XsCOofJ 

(26) 


ton 


As  before,  tlu'  doubh'  series  in  this  equation  reduces  to  a 
single  series  if  lb(>  discussion  is  limited  to  tlie  intorfen'iice  at 
the  center  lim'  of  the  tunnel.  For  points  on  the  center  line, 

Jn(X,,.V)  =  Jn(0)  =  ] 

(X,,i/)  =  J,„(0)  =  0  form^  1 
and  tlie  streamwise  indticed  velocitv  becomes 


/ _ r  dr]  I 

''  "  27r/'-^^ 


sin  0nt/(/(XscO(,) 


Jid(Xsr) 

From  the  knowti  relations  foi'  the  Bessid  functions 

</i/(X,,t0n)  =  — t7|(X,,a)n) 

so  that  equation  (27)  tnay  be  written 

/_  _t  t/77 sin  0(1  ^7, (X.s'ton) 
~‘2Trr-^,  Jn-(X,,r) 


(27) 


(28) 


(20) 


As  required  by  etiuation  (4),  the  summation  with  resiject  to 
.xf  in  this  equation  e.xtends  over  all  positive  roots  of  the 
equation 

'A,  l7|(X,v/')  — d  (50) 

As  the  ne.\t  step,  consider  a  pair  of  symmetrically  placed 
i'lementary  vortices  composed  of  a  vortex  of  circulation  —  T' 
at  tlu‘  point  oio,  8,,  and  a  vortex  of  circulation  +  F'  at  the 
jjoint  CO, I,  0;^  Fi-om  ('quation  (29),  tlu'  streamwise  velocity 
induci'd  at  a  point  on  the  center  line  of  the  tunnel  by  this 
vortex  jiair  and  the  accompanying  trailing  vortices  is 


T'dr]  ■ 


X  -  - 


sin  0„,/i(X,,coo) 

which  may  also  be  written 

2t''(0|)Sili  0„  c/r)  ^  e'''V,(XsCo„) 


(51) 


(52) 


The  expression  (2r'con  sin  9,,)  which  ap])ears  in  this  cqtiation 
is  the  product  of  the  vortex  strength  and  the  distance  bctwc'cn 
the  vortices. 

Now  let  tin'  distance  betweem  tlu'  vortices  tend  to  /-('ro 
while  the  vortex  sti'cngtb  increases  in  such  a  way  that  the 
product  (2r'coo  sin  0o)  redains  a  constant  value  ji.  Tin' 
residt  in  the  limit  is  an  elementary  vortex  doublet  of  strcuigth 
fi  and  span  dg  at  the  point  000=17  on  the  y-axis.  Tin' 
streamwise  velocity  induced  on  the  center  line  of  the  tunmd 
by  this  elementary  spanwisc  doublet  and  the  accompanying 
trailing  vortex  doublets  is  then 


(U  c/17  e>‘»L/,(X,i7) 

27rr-i7  Jo‘(X,,/') 


(55) 


As  before,  the  infinite  series  in  this  equation  is  rapidly 
convergent  for  large  lu'gatiAU'  values  of  x,  but  tlu'  convergence 
is  slow  for  small  negative  value's  and  is  nonexistc'ut  wlu'ii 
x=0.  Once  sigain,  bowevc'r,  the  series  can  lx-  expressc'd  as 
a  combination  of  elementary  functions  and  a  power  series 
which  is  readily  applied  to  the  probb'in  at  band.  The 
details  of  the'  transformation  are  given  in  Appc'iidix  A. 
By  meiins  of  the  final  result,  ('(|uatioJi  (.53)  can  be  writtc'ii 

..  /  M  dl7  r  (  1  )^y2(<+;i-M)17^br^'’  "1 

~  /-!(7'+ 1  )!(2il)!2A‘+q,2c-  +  2;M  ,  J 


(54) 


f"  /X.. 


where  the  double  summation  extends  over  idl  integral  values 
of  7' and  p  from  zero  to  infinity.  The  coefficient  y2(*:+?,+  i)  = 
Hif  is  given  by  the  integral 

fdh 

'(2/+1)irj„  Idif) 

'I'lie  numerical  values  of  this  integral  for  /=!,  2,  5,  4  are 
evaluated  in  Appendix  A. 

'The  induced  velocity  for  a  doublet  spanning  the  tunm'l  is 
now  readily  found  by  taking  the  doublet  strength  ,«  constant 
across  the  span  and  integrating  equation  (54)  with  respect 
to  17  from  — /■  to  +/•.  This  gives  finally 


/ _ M  r  ^  l)^M2(A:  f;M  I)T 

■'  2,rFLx;Xd  /-  f^,^»k\[k+m2v)\{2k+\)2 


1)2-V'’_ 


(56) 


In  the  integration  across  the  tunnel,  idl  the  trailing  vortices, 
of  course,  disappear. 

It  is  apparent  from  the  symmetry  of  the  jji'oblem,  that  the 
streamwise  velocity  induced  by  a  doublet  spanning  the  tun¬ 
nel  must  be  an  even  function  of  the  variable  i.  Equation 
(56),  which  was  derivxed  for  negative  values  of  x,  is  seen  to  lie 
such  a  function  and  is  thus  applicable  to  posithm  values  of 
the  variable  as  well. 

The  values  of  c/  for  vanishingly  small  vahu's  of  x,  that  is, 
at  the  position  of  the  doublet,  is  then  found  from  ecpiation 
(56)  by  expanding  the  first  term  in  ascending  powers  of  xjr 
and  discarding  all  terms  containing  second  powers  and 
higher  and  by  retaiinng  only  the  p  =  Q  terms  of  the  double 
series.  This  gives. 
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After  substitution  of  the  numoTical  values  for  the  coefficients 
M2(i+i)  from  equations  (A33)  of  Appendix  A,  this  equation 
becomes  to  an  accuracy  of  three  significant  figures 


/I  I  f. 3 56/2 
2irJ-  2ivr^ 


f38) 


The  first  term  of  equation  (38)  is  the  vcdocity  induced  by  a 
doubhd  of  infinite  span  in  a  field  of  unlimited  extent.  The 
remaining  t('i’m  thei-efoi'c  ii'presents  the  (‘fl'ect  of  interference 
between  flie  doubhd  and  the  tuinud  wall.  Thus  the  net 
result  of  tlu'  interference  between  tlie  airfoil  thickness  and 
the  tunnel  wall  foi'  the  incompressible  fluid  is  to  increase  the 
effective'  stream  velocity  at  the  position  of  the  airfoil  by  the 
amount 


A,F' 


.356/2 

27r/" 


(39) 


In  any  j/articular  case,  /2  is  again  equal  to  the  strength  of  the 
doublet  used  to  re]5resent  the  given  airfoil. 

The  result  of  equation  (39)  can  be  modified  for  the  eff/'ct 
of  fluid  compressibility  by  the  method  of  reference  4.  In 
this  case,  the  modification  is  most  conveniently  performed 
by  means  of  Method  1  (reference  4,  pp.  3-5)  which  compares 
the  comjiressihle  and  incompressihle  flows  for  a  given  airfoil 
of  unaltered  shape'  and  size.  By  this  method,  it  is  readily 
shown  that  the  streainwise  velocity  induced  in  the  incom¬ 
pressible  fluid  at  any  point  on  the  center  line  of  a  tunnel  of 
radius  r  is  times  the  cori'esponding  velocity  at 

the  saiiK'  point  in  an  incompressible  fluid  in  a  tunnel  of  radius 
r-y/l  —  (M')-.  Hei'e  M'  is,  as  before,  the  Mach  number  in  the 
undisturbed  stream.  The  increment  in  axial  velocity  in  the 
compri'ssible  case  is  thus 


AiV’ 
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(40) 


Comparison  of  etpiations  (24)  and  (40)  shows  that,  if  no 
wake  is  prc'sent,  a  symmetrical  aiifoil  spanning  a  closed- 
throat  circular  tunnel  of  radius  r  experiences  at  its  midspan 
section  the  same  increase  in  axial  velocity  as  would  bc'  ex¬ 
perienced  by  the  same  airfoil  in  a  closed-throat  two-dimen¬ 
sional-flow  tunnel  of  height 


Wake  effect. — It  is  shown  in  general  terms  in  reference  5 
that  the  interference  between  the  wake  of  a  body  and  the 
walls  of  a  closed-throat  wind  tunnel  gives  rise  at  the  position 
of  the  body  to  a  velocity  increment  and  a  streamwise  pres¬ 
sure  gradient  which  are  not  present  in  free  air.  This  is  true 
for  any  type  of  body  and  any  shape  of  tunnel  test  section. 
The  magnitude  of  this  velocity  increment  and  pressure 
gradient  in  the  case  of  an  airfoil  spanning  a  closed-throat 
rectangular  tunnel  can  be  determined  approximately  by 
replacing  the  wake  hy  the  flow  from  a  suitable  fluid  source 
and  the  tunnel  walls  by  an  infinite  system  of  image  sources. 
In  the  case  of  the  airfoil  spanning  a  closed-throat  circular 
tunnel,  this  treatment  is  no  longer  possible,  since  no  system 
of  image  sources  is  known  which  will  satisfy  the  boundary 
conditions  at  the  tunnel  wall.  .V  more  comple.x  method  of 
analysis  could  conceivahly  be  devised  for  this  case;  however, 
since  the  calculation  is  highly  approximate  even  in  the  case 
of  two-dimensional  flow,  such  an  analysis  does  not  appear 
warranted.  For  present  purposes  it  is  probabh^  sufficient  to 
assume  that  the  midspan  section  of  the  airfoil  in  the  circidar 
tunnel  experiences  the  same  velocity  increment  and  pressure 
gi-adient  as  a  result  of  the  wake  interference  as  does  the  same 
airfoil  in  a  rectangular  tunnel  of  a  height  defined  by  equa¬ 
tion  (41).  This  assumption  leads  to  the  simplest  expression 
for  the  final  correction  to  the  measured  drag  coefficient  and 
should  give  results  which  are  reasonably  accurate.  If  it  is 
assumed  that  the  center  of  the  wake  lies  in  a  horizontal  plane 
containing  the  diameter  of  the  tunnel,  it  follows  from  con¬ 
siderations  of  symmetry  that  the  wake  interference  does  not 
contribute  to  the  vertical  induced  velocity  a/  at  the  airfoil. 

It  has  already  been  indicati'd  that  the  interference  as¬ 
sociated  with  the  camber  of  the  airfoil  has  no  ('fl'ect  upon  the 
stream  velocity  at  tlic  model.  The  total  increase  in  velocity 
for  the  complete  aiifoil  in  the  circular  tunnel  is  thus  given  by 
the  sum  of  the  increments  caused  by  the  thickness  and  tlie 
wake  of  the  aii'foil.  In  reference'  5  it  is  shown  that  for  tlu; 
analagous  case  of  tlu'  airfoil  in  the  rectangidar  tunnel,  the 
true  velocity  1’  at  the  position  of  tlu'  airfoil  may  finally  l/e 
written 


(42) 
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.y  r  =  1.5.58r 


or,  in  tc'i'ins  of  the  tuniu’l  diametc'r, 

h,  =  0.779d  (41) 

The  foregoing  l  esult  greatly  simplifu's  the  determination  of 
the  tnu'  strc'am  conditions  at  the  position  of  tlu'  airfoil  in  the 
circidar  tunnel,  since  tlu'  lu'ct'ssary  (’quations  foi’  the  ri'c- 
tangular  tunnel  are  alr('ady  known. 

Consideration  of  the  symmetry  of  the  system  fornu'd  by  a 
base'  pi'ofile  spanning  the  middle  of  a  circular  tunnel  indicatc's 
that  the  int ('reference  Ix'tvveen  the  wall  and  the  airfoil  thick- 
lU'ss  does  not  infliu'nce  the  vertu'al  induced  velocity  v^'  at 
any  point  on  the  aii'foil.  Similarly,  the  airfoil  thickness  has 
no  ('ffe('t  upon  the  streamwise  pressure  gradi('nt  in  the  tunnel 
at  th('  position  of  tlu'  airfoil. 


where  u  and  t  are  factors  (h'pendent  upon  tlu'  siz('  of  the  aii  - 
foil  relative  to  the  tunnel,  A  is  a  factor  dep('n(lent  upon  the 
shape  of  the  base  profile,  and  c/  is  tlu'  drag  co('ffici('nt  of  the 
airfoil  as  measured  in  the  tunm'l.  Tlu'  first  correction  term 
in  this  equation  repr('sents  the  vi'locity  increment  caused  by 
the  airfoil  thickness  and  is  found  by  substituting  the  proper 
value  for  the  equivah'ut  doublet  str('ngth  in  ('(pialion  (24). 
Thes('con(l  cori-(X'.tion  U'rm  r('pr('S('nts  the  v('locity  incr('ment 
associated  with  tlu'  wake  of  the  airfoil. 

The  factors  a  and  t  in  equation  (42)  are  defined  by 
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where  (c/A)  is  the  ratio  of  the  airfoil  chord  to  the  tunnel 
height.  ;\n  analytic  expression  for  A  is  given  in  equation  (3) 
of  reference  5.  Values  of  A  for  a  number  of  base  profiles  are 
given  in  table  I,  which  is  reproduced  from  this  reference. 

If  it  is  assumed  that  the  height  of  the  eqinvalent  rectan- 
gulai'  tunnel  with  regard  to  the  wake  interference  is  the  sanie 
as  that  given  by  equation  (41)  for  the  thickness  interference, 
the  true  velocity  in  the  ciretdar  tunnel  is  found  simply  by 
sul)st  ituting  h  ,  from  equation  (41)  for  h  in  th('  factors  a  and  t 
of  ('quation  (42).  The  true  velocity  at  the  niids])an  R(‘c1ion 
of  ail  airfoil  spanning  a  circular  tunnel  is  thus 


V- 


3/2  A  0-2  + 


l+0.4(d/')^'  ,} 


(43) 


a.nd  0-2  are  defini'd  by 

,,-0.301  (;,) 

(40) 

.,™0.339(')' 

(17) 

.V  correction  to  the  stream  velocity  inqilies  corrections  also 
to  th(‘  stream  dynamic  pressunu  Reynolds  nundx'i-  and 
Mach  number.  These  corrections  for  an  airfoil  spanning  a 
rectangidar  tunnel  have  been  determined  in  reh'i-enci*  5  on 
tJu'  basis  of  the  assumption  that  the  flow  is  adiabatic 
d’tie  eoiTcsponding  corrections  for  the  cii'cular  tunnel  can 
be  found  by  replacing  tlu'  factors  t  and  cr  in  equations  (29), 
(32),  and  (33)  of  reference  5  by  the  factors  and  cr>  of  the 
prc'seiit  paper.  The  true  dynamic  pi'essure  q,  Rcwnolds 
number  R,  and  Mach  number  Af  at  tlie  niidspan  section  in 
the.  circular  tunnel  are  thus  related  to  the  corresiionding 
(juanlities  in  the  undisturbed  stream  (denoti'd  by  jirinies)  liy 
t he  equations 


<1-  <1' ,  I  + 


(  '  ^  [1--(A/')T'' 


A  (T.J  H  - 


(■IS) 


R--/r\i  + 


-O.liAf'Y  [l-l).7(A/')-’][l -I  (),4(.l/')d  d 

-(iU')T-  '  ]-{M'Y  \ 

149) 


I  '  0.2(.l/'')-  [1  -L0.2(;\/')-][l  \  0A(M'Y]  ,] 

[l-(d/')T''-  l-(d/')"'  -  ’-T./  ^ 

(AO) 


Numiu’ical  values  of  the  functions  of  A/'  wliich  a])i)ear  in 
these  equations  are  given  in  tabh'  11,  which  is  i-e]n'oduced 
from  reference  5. 

At  low  Mach  numbi'rs,  the  tiM’ius  containing  t-,c/  in  tlu' 
c(|uations  for  the  corrected  sti'i'am  cbai’aeteristics  are  usually 
negligible  as  couqrarc'd  with  the  terms  containing  Acri.  At 
high  Mach  numbi'rs,  howevi'r,  where  the  ding  coeflicienl  is 
very  large,  the  ti'rms  with  roC,/  ])redominate. 

RKLATIONS  BKTWUKN  AIRFOIL  CH AKACTKKI.STICS  IN  TUNNUL 
AND  IN  FREE  AIR 

d'hi'  characteristics  of  the  airfoil  in  fri'e  air  are  now  readily 
determined  in  ti'rms  of  the  charactcj'istics  at  the  midspa.n 


section  in  the  tunnel.  It  is  simply  necessary  to  apply  the 
results  of  the  preceding  sections  to  the  relations  already 
derived  in  reference  5  for  the  airfoil  spanning  a  rectangular 
tunnel. 

Briefly,  the  method  of  reference  5  relates  the  section  char¬ 
acteristics  in  the  tunnel  at  an  undisturbed  stream  velocity 
V'  to  the  characteristics  in  an  unconfined  stream  having  a 
velocity  equal  to  the  true  velocity  V  which  exists  at  the 
position  of  the  airfoil  in  the  tunnel.  The  relation  is  ob- 
lained  on  the  basis  of  equal  values  of  the  so-called  cotangent 
component  of  lift  in  the  tunnel  and  in  free  air,  this  being 
necessary  to  assure  that  the  essential  character  of  the  ])r('s- 
sui'c  distribution  over  the  airfoil  is  the  same  in  both  cases. 
By  this  ])rocedure  corrections  are  derived  which  may  be  ap¬ 
plied  to  simultaneously  measured  lift,  moment,  and  drag  co¬ 
efficients  and  angle  of  attack  in  the  tunnel  to  obtain  the  cor¬ 
responding  quantities  in  free  air.  These  corrections  appear 
as  functions  of  the  factors  A  and  a,  of  the  produet  rc/ ,  and 
of  the  Mach  number  M'  of  the  undisturbed  stream.  The 
correction  to  the  angle  of  attack,  wdiich  arises  out  of  tlu' 
interference  eflects  associated  with  camber,  is  proportional 
to  cr  and  independent  of  A  and  tc/.  The  correction  equa¬ 
tions  for  the  lift  and  moment  coellicients  contain  correspond¬ 
ing  terms  proportional  to  a  alone,  together  with  terms  which 
depend  upon  the  thickness  and  the  wake  effects  and  are  pro- 
jjortional  to  the  pi'odncts  Ac  and  tc/.  The  correction  to  the 
drag  coefficient  appears  as  two  terms,  proportional  to  Aa 
and  tc/,  respectively.  The  term  proportional  to  Ao-  is  in 
this  case  composed  basically  of  two  parts,  one  due  to  the 
thickness  effect  and  one  due  to  the  w'ake  effect. 

The  coi'ri'ction  eejuations  foi'  thi'  airfoil  spanning  a  circular 
tunnel  can  be  di'rivi'd  diri'ctly  by  modifying  tlie  ('{piafions  of 
I'eh'rence  5  in  accordance  wdth  the  results  of  the  preceding 
sections.  Since  the  terms  containing  <r  exclusive  of  A 
appear  as  a  result  of  the  camber  effect,  the  tunnel  height  h, 
in  such  terms  must  be  replaced  by  0.843d  as  required  by 
e(|uation  (23).  In  the  terms  which  depend  upon  the  thick¬ 
ness  and  wake  effects  and  ari'  distinguished  by  the  products 
A<7  and  tc/ ,  tlu'  quantity  h  is  replaced  by  0.779d  in  accord¬ 
ance  with  equation  (41).  This  involves  the  assumption 
already  mentioned  that  the  height  of  the  I'quivalent  rec¬ 
tangular  tunnel  with  ri'ganl  to  tlu'  wake  effect  is  the  same 
as  that  calcidated  for  the  thickness  effect. 

As  in  reference  5,  the  free-air  lift,  quarter-cbord-moment, 
and  drag  coefficients  referred  to  the  true'  dynamic  jiri'ssuri'  q 
are  di'iiotcd  by  ihe  conventional  sjunbols.  Tlu'  corri'spond- 
ing  quantities  nu'asurcd  in  tlu'  tunnel  and  referred  to  the 
apparent  dynamic  pressure  q'  are  denoted  by  the  same 
symbols  with  jn'imes  adiled.  The  final  equations  for  the 
corrected  aerodynamic  coefficients  are  then 

<^1  2-{M'Y  . 

y  i-(Af')^  [r-(A/')^P  ' 

[2-(A7')^][l+0.4(AT)^] 

1-(A/')-  ^  ) 


c 


2— _ 

[2-^(Ar)^][l+0.4(A/')1 

~  ]~(A1'Y  """  r  4[1-(A/0^1 


(r,2) 
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0-0^1  [l_(iJ4-')2]3/2  Aff2 

[2-(M')-][l+0.4(M')^]  ,) 

iiiid  (ho  oon-os|)oti(ling  anglo  of  attack  in  dogriM's  is 


wluvfo  till'  facl-oi'  (T|  is  giv(‘n  liy 

,,  =  (1^89  (0‘ 


(5:5) 


(r,4) 


(55) 


and  (ho  factors  to  and  tro  ari'  as  alrc'ady  defined  in  equations 
(40)  and  (47).  Numerical  values  of  the  compressibility 
factors  which  appeal'  in  these  equations  are  given  in  table  IT. 
The  cori'ected  quantitii's  correspond  to  the  true  Reynohls 
number  and  true  Mach  number  as  given  by  equations  (41)) 
and  (50). 

From  a  rigorous  standpoint,  the  foregoing  corrections 
aiipiy  only  to  data  obtained  from  chordwise  pressure  dis¬ 
tributions  at  the  midspan  section  of  tlio  airfoil.  Actually, 
as  baa  ah'(>ady  been  pointed  out  in  the  discussion  of  camber 
effect,  the  exjierimental  chordwise  pressure  distribution  at 
any  given  angle  of  attack  is  sensibly  constant  across  the  span. 
Tb(!  corrections  should  therefore  be  applicable  with  sufficient 
accuracy  to  data  obtained  fi'oni  pressure  distributions  at  any 
spanwise  station. 

K('f('rence  5  also  includes  a  method  for  cori'octiiig  experi- 
nuuital  cliordwise  ju'C'ssurc  distributions  to  free-air  condi¬ 
tions  in  the  case  of  an  airfoil  spanning  a  rectangular  tunnel, 
d'he-  same  procedure  may  be  applied  to  pressure  distributions 
ovi'i'  an  airfoil  spanning  a  circular  tunnel  if  the  factor  t 
is  re]daced  by  to  and  the  factor  <r  by  a,  wlierever  it  appi'ars 
alone  and  liy  ff2  wher(‘  it  appears  in  the  pi'oduci  Act. 

CHOKING  AT  HIGH  SPEEDS 

As  explained  in  reference  5,  for  tests  of  a  model  in  any 
closed-llii'oat  wind  tunnel,  then'  is  some  value  of  (be  Macb 
number  A/'  of  the  undisturbed  stream  which  cannot  be 
exet'eded  irres]iective  of  the  powei'  input  to  the  tunnel. 
This  follows  from  the  fact  that  at  high  speeds  the  combina¬ 
tion  of  model  and  wind  tunnel  acts  essentially  as  a  converg¬ 
ing-diverging  nozzle,  and  the  flow  in  the  tnnnel  exhibits  tlie 
characteristics  of  the  flow  in  such  a  nozzle.  'Plius,  at  some 
.Mach  numbi'r  less  than  unity  in  tlu'  undisturbed  stream, 
sonic  velocity  is  attained  at  all  jioints  across  a  setUion  of  the 
tunnel,  usually  in  the  vicinity  of  the  model.  When  this 
occurs,  increased  power  input  to  the  tnnnel  serves  merely 
to  extoiid  the  region  of  su]jersonic  flow  downstream  of  this 
sonic  section  and  has  no  ('fi'ect  iqion  tlu'  velocity  of  the  stream 
ahead  of  tbi!  model.  The  tunnel  is  then  said  to  be  “choked”, 
and  the  Mach  number  AI'  of  the  undisturbed  flow  ahead  of 
I  lie  aii'foil  has  its  maximum  attainable  value,  d'his  valtie  is 
described  as  tlu'  appai'cnt  choking  .Mach  numbi'r,  tin',  word 
“ap]uu'en(”  being  used  to  ditferentiate  this  value'  from  the 
coi'i'esponding  fri'c-air  Mach  number  A/  which  would  be 
computed  from  equation  (50). 


If  it  is  assumed  that  the  section  of  sonic  velocity  is  coin¬ 
cident  with  the  section  of  minimum  area  between  the  model 
and  the  tunnel  walls,  the  apparent  choking  Mach  number 
can  be  obtained  from  elementary  considerations  of  uni- 
dimensional  adiabatic  flow,  as  shown  in  reference  5.  Foi' 
the  present  case  of  a  constant-chord  airfoil  spanning  a  cir¬ 
cular  tunnel,  the  apparent  choking  Mach  number  A/T*  if 
finally  defined  for  air  (7=1,4)  by  the  relation 


4 

TT 


1  + 


(MT;d- 

6 


■7 


(50) 


where  4  is  the  “effective”  thickness  of  the  airfoil  and  d  is,  as 
hefore,  the  diameter  of  the  tunnel.  A  graph  of  this  relation 
is  given  in  figure  4.  As  a  matter  of  interest,  the  results  are 
shown  for  the  super.sonic  as  well  as  the  subsonic  flow  regime, 
d'he  region  above  the  curve  represents  an  impossible  state 
of  flow. 

In  estimating  the  apparent  choking  Mach  number  in  any 
practical  case  it  is  necessary  to  replace  the  effective  thickness 
ig  by  the  projected  thickness  tp  of  the  airfoil  normal  to  the 
direction  of  flow.  As  indicated  in  reference  5,  this  procedure 
leads,  in  the  case  of  the  subsonic  wind  tunnel,  to  an  over¬ 
estimation  of  M'cft  because  it  neglects  the  possible  con¬ 
traction  of  a  portion  of  the  stream  aft  of  the  airfoil  as  well 
as  the  elfect  of  the  airfoil  boundary  layer. 


Choking  Mach  Numbenj  M'ch 


FKjritK  1.-  Ohokiiig  Mtich  nunibcr  :ts  (Icfi'iinini'tl  by  thickness. 
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The  importance  of  the  boundary  layer  and  the  accom¬ 
panying  drag  with  regard  to  tunnel  choking  is  pointed  out 
in  reference  5,  where  the  apparent  choking  Mach  number  is 
calculated  for  a  flat  plate  at  zero  angle  of  attack  in  a  two- 
dimensional-flow  wind  tunnel.  Since  the  projected  thickness 
for  the  plate  is  zero,  the  unidimensional  theory  would  indicate 
that  no  choking  occurs.  Actually,  because  of  the  fact  that 
the  ])late  experiences  drag,  choking  iloes  take  place.  Similar 
considerations  hold,  of  course,  for  a  flat  plate  spanning  a 
cii'cnlar  tunnel,  in  this  case  tlu'  apparent  choking  Mach 
number  foi’  aii'  (7=1.4)  is  given  by  the  equation 


no  way  is  known  to  combine  the  thickness  and  drag  eff(‘cts 
in  a  single  calculation  as  should  logically  be  done. 

It  should  be  noted,  as  pointed  out  in  detail  in  reference  5, 
that  the  flow  in  a  tunnel  at  choking  does  not  correspond  to 
any  flow  in  free  air.  Furthermore,  for  a  range  of  Mach 
numbers  just  below  choking,  where  the  flow  is  influenced  to 
any  extent  by  the  restrictions  which  finally  promote  choking, 
any  wall-interference  correction  is  of  doubtful  accuracy. 
This  is  particularly  true  if  the  model  is  at  an  appreciable 
ai>gle  of  attack  so  that  sonic  velocity  is  attained  across  the 
stream  on  one  side  of  the  airfoil  before  it  is  on  the  other. 


1  +  / 

2.8(M'„y  I  ^  V 


f-(A/h„)^  "IM 
J  +  1.4(Af',,)^J  i 


(57) 


A  graph  of  this  relation  is  given  in  figuii'  5.  Tin'  elh'ct  of 
drag  on  choking  for  supersonic  as  well  as  subsonic  wind 
tunnels  is  shown.  It  can  be  dcunonstrated  that  the  points 
on  the  curve  coi'respond  to  a  Alach  number  of  unity  in  tin* 
flow  far  downstream  of  the  model  whei'e  the  wake  has  spread 
completely  to  the  tunnel  wall.  Points  above  the  curve 
repi'csent  impossible  conditions  of  flow.  In  most  cases 
(uicountered  in  subsonic  tunnels,  the  apparent  choking 
Mach  number  determined  by  the  thickiu'ss  of  the  airfoil  and 
defined  by  equation  (56)  is  usually  the  lower.  For  very  thin 
airfoils  at  small  angles  of  attack,  howevei',  the  valiu'  of  M'c), 
given  by  equation  (57)  can  have  the  lower  value.  At  present 


Choking  Mach  Number,  M'ct, 


EXPERIMENT 

The  experimental  investigation  was  initiatc'd  for  two  r('a- 
sons:  (1)  to  determine  the  spanwise  distribution  of  lift  over 
an  airfoil  spanning  a  closed-throat  circular  tuniud,  and  (2) 
to  examiiu'  the  validity  of  the  theoretical  intcrh-rence  (nr- 
rections  deriv<‘d  in  the  preceding  analysis.  As  has  bc-en  pre¬ 
viously  mentioned,  the  development  of  the  theoi'etical  rela¬ 
tions  requires  a  knowh'dge  of  the  variation  in  lift  over  the 
span  of  the  airfoil.  Since  no  theoretical  or  experimental 
evidence  regarding  this  matter  was  available,  the  spanwisc' 
variation  in  lift  was  investigated  experimentally  for  an 
NACA  4412  airfoil  for  two  ratios  of  airfoil  chord  to  tunnel 
diameter.  The  results  of  those  tests  are  also  directly  appli¬ 
cable  to  the  examination  of  the  validity  of  the  theoretical 
coiTcction  equations. 

The  experimental  woi'k  was  peilormed  in  a  low-turbulence, 
nonreturn-type  wind  tunnel  w  ith  interchangeable  throat  sec¬ 
tions  of  14-  and  8-inch  diameter.  The  two  choid-diameter 
ratios  were  obtained  by  testing  the  same  airfoil  in  each  throat 
section.  Since  the  airspeed  was  held  constant  throughout 
the  tests,  this  arrangement  permitted  the  Reynolds  number 
and  the  Mach  number  to  be  duplicated  simultaneously  for 
the  two  chord-diameter  ratios.  In  this  manner  the  eflects 
of  any  variation  in  these  parameters  were  eliminated  from 
the  tests. 

The  NACA  4412  airfoil  was  used  because  a  model  of  suit¬ 
able  size  was  already  available  ideally  equipped  for  pressure- 
distribution  tests.  The  model,  which  is  desci-ibed  in  refer¬ 
ence  16,  was  of  5-inch  chord  and  30-inch  span.  This  chord, 
together  with  the  two  throat  diameters,  gave  chord -diameter 
ratios  of  0.357  and  0.625.  In  the  tests,  the  airfoil  extended 
through  the  walls  of  the  tunnel  and  was  clamped  in  tight- 
fitting  support  blocks  which  prevented  any  leakage  of  air  at 
the  walls.  The  54  pressure  orifices  located  around  the  sur¬ 
face  of  the  midspan  section  of  the  model  w^ere  connected  to 
a  multiple-tube  manometer  for  measurement  of  the  pressure 
distribution  over  the  airfoil.  To  secure  as  accurate  pressure- 
distribution  data  as  possible,  alcohol  was  used  as  the  manom¬ 
eter  fluid  and  the  liquid  heights  w'ere  recorded  photographi¬ 
cally. 

Pressure-distribution  records  were  secured  at  each  of  eight 
angles  of  attack  from  —4°  to  15°  at  a  Reynolds  number  of 
approximately  450,000  and  a  Maidi  number  of  approximately 
0.2  with  the  model  mounted  in  both  the  14-inch  and  the  8- 
inch  diameter  throats.  The  spanwdse  distribution  of  lift  w  as 
determined  for  each  angle  of  attack  by  sliding  the  pressure 
orifices  laterally  from  one  wall  to  the  otlu'r  and  rc'cording  tlu' 


Fioure  .h— Choking  Mach  nuni))er  ay  determined  by  airfidl  drag. 
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indicated  pressure  distributions  at  a  number  of  spanwise. 
stations.  The  cbordwise  pressure  distributions  were  plotted 
aiul  meebanically  integrated  to  obtain  lift  and  quarter-chord 
moment  coefficients.  No  drag  coefficients  were  obtained 
because  the  experimental  installation  did  not  permit  balance 
nu'asurements  to  be  made  and  wake  surveys  were  not  feasible. 

By  testing  the  airfoil  in  both  erect  and  inverted  attitudes 
the  inclination  of  the  air  stream  with  respi'ct  to  the  tunnel 
axis  was  determined  for  each  throat  section.  The  stream 
angle  was  found  to  be  +0.45°  for  tlu'  14-inch  throat  and  0° 
for  tlu'  8-inch  throat,  (hu-rections  have  Ixaui  applied  to  all 
angle's  of  attack  for  the  nu'asurc'd  angularity. 

The  sj)anwis('  distribution  of  lift  coefficie'iit  uncorrectcil  for 
tunnel-wall  inte'i'ference  is  shown  for  the  two  chord-diameter 
ratios  in  figure's  6  anel  7  in  which  lift  e'oefficients  at  x^arious 
angles  of  attai'k  are  plotted  as  a  function  of  the  spanwise 
location  of  the  me'asurenient  plane. 

Curves  of  lift  coefficie'nt  against  angle  of  attack  for  the  twe) 
e'hord-elianie'ter  ratios  are'  shown  uncorrecte'd  for  tunne'l-wall 
inte'i'fcrene'e'  in  figure  8  (a).  The'  results  given  pe'rtain  to  the 
se'ction  of  the'  aii'foil  at  the'  center  line  of  the  tunnel.  The 
e'orres]ioneling  e'urve's  cori-e'cted  for  wall  interference  by 
nu'ans  of  e'epiations  (51)  anel  (54)  are  shown  in  figure  8  (b). 
In  applying  tbe  corre'ctious,  the'  term  containing  rCa  was 
ne'e'e'ssarily  omitteel  as  no  measure'uie'nts  of  drag  we're'  made'. 
Tor  the'  value's  of  Ca  to  be'  e'xpecte'el  in  such  te'sts,  howe've'r, 


Fic.ure  G.— Spiinwise  lift  distiibutioii  for  lhi“  N  A(\\  1112  airfoil.  c/d=.357. 


Dis  fence  from  the  funnel  cenferline  in  percenf  radius 
Fhiuhk  7.— Spanwise  lift  (ILstribution  for  the  XAOA  4412  airfoil,  c/(/=.G2rL 


this  te'rm  woulel  be'  lu'gligible  in  comparison  with  the'  re'main- 
ing  terms  so  that  this  omission  does  not  affee't  the  final  results. 
.For  purpose's  of  e'omparison,  section  lift  e'harae'teristieis  as 
olitained  by  Pinkerton  from  tests  of  a  finite'-span  re>ctanguhir 
airfoil  in  the  Langh'y  variable-elensity  wind  tunned  (reference' 
17)  are  also  shown.  These  data  correspond  to  an  effective 
Re^ynolels  number  of  450,000  anel  are  thus  directly  comparable 
to  the'  I'e'sults  of  the  jiresent  test. 

In  figure'  9  (a)  curve's  of  quarter-e'horel  mome'iit  e'ot'ffie'ie'iit 
against  lift  e'oe'fhcient  are  shown  une'orree'te'el  for  tnnnel-wall 
intcrfei'cnce'  for  both  chorel-diamcter  ratios.  The  same'  data 
are  plotteel  in  figure  9  (b)  after  e'orrection  for  wall  intei'fe'rene'e 
by  means  of  eepiations  (51)  and  (52).  Also  shown  for  com- 
jiarison  are  the  corresponding  data  from  reference  17. 

DISCUSSION 

An  examination  of  figures  6  and  7  reveals  the  previously 
mentioned  fact  that  there  is  no  appreciable  variation  in  lift 
over  the  span  of  the  airfoil  at  all  angles  of  attack  up  to  those 
closely  approaching  the  stalling  angle.  This  observation 
holds  for  both  chord-diameter  ratios.  In  the  vicinity  of  the 
stall  a  spanwise  variation  in  lift  appears  which  becomes  pro¬ 
gressively  more  erratic  as  the  angle  of  attack  is  increased. 
As  might  be  expected,  this  variation  becomes  apparent  at  a 
lower  angle  in  the  case  of  the  larger  chord-diameter  ratio. 
The  rcsidts  of  figures  6  and  7  corroborate  the  conclusion  of 


Seciion  piichmg  momeni  coeffidenf^  C^»  Section  /iff  coefficient, 
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Appendix  B  for  the  particular  case  of  the  airfoil  spanning  a 
circular  tunnel. 

From  figure  8-  (b),  it  is  seen  that  the  corrected  lift  curves 
for  the  two  chord-diameter  ratios  agree  almost  exactly  with 
one  another  except  at  angles  near  the  stall.  Below  the 
vicinity  of  the  stall  the  corrected  data  coincide  with  the 
results  of  reference  17  except  for  a  constant  angular  displace¬ 
ment  of  approximately  0.2°.  In  reference  16,  Pinkerton 
estimates  that  his  values  for  the  angles  of  attack  may  he 
too  large  by  a  constant  error  of  approximately  0.25°  because 
of  a  possible  eiTor  in  the  assumed  direction  of  the  stream. 
It  is  thought  that  the  angles  of  attack  of  the  present  experi- 
mer.ts  are  accurate  to  within  ±0.1°.  These  limits  of 
accuracy  arc  sufficient  to  account  completely  for  the  apparent 
angular  displacement. 

In  the  region  of  the  stall,  the  corrected  lift  curves  for  the 
two  chord -diameter  ratios  do  not  mutuallv  coincide,  hut 
the  data  for  the  chord-diameter  ratio  of  0.357  agree  witli 
Pinkerton’s  results  within  2  percent.  As  previously  men¬ 
tioned,  Pinkerton’s  tests  were  made  with  a  finite-span 
r('ctangular  airfoil,  for  which  the  cross-span  variation  in 
lift  is  necessarily  large.  It  is  not  to  be  expected  that  the 
determination  of  maximum  section  lift  from  such  tests 
would  be  as  accurate  as  from  tests  of  a  through  model,  foi' 
which  the  cross-span  lift  variation  is  small. 

It  is  seen  from  figure  9(b)  that  the  corrected  moment 
curves  agree  satisfactorily  with  each  other  and  with  the  re¬ 
sults  of  reference  17. 

In  summary,  for  angles  of  attack  below  those  in  tlu*  I'cgion 
of  maximum  lift,  the  rissults  presented  in  figui'os  8  and  9 
demonstrate  the  validity  of  the  theoretical  lift,  moment,  and 
anglo-of-attack  corrections  foi'  low  Alach  numbers  and  cliord- 
diametcu'  ratios  up  to  at  least  0.625.  For  angles  in  tlu'  vicin¬ 
ity  of  maximum  lift,  the  corrections  are  not  strictly  ajiplica- 
ble  u])  to  such  a  large  chord-diameter  ratio.  The  j'csults  of 
the  present  ti'st  indicate  that  an  accurate  determination  of 
maximum  lift  can  be  made  with  a  chord-diameter  ratio  at 
least  as  high  as  0.35.  An  evaluation  of  the  accuracy  of  the 
correction  equations  at  high  Mach  numbers  is  not  possible 
on  the  basis  of  the  experimental  evidence  availabh'  at  pi'esent. 
It  is  to  be  expected,  however,  that  the  maximum  p(>rmissible 
chord-diameter  ratios  will  decrease  as  the  Mach  iiumbcu-  in¬ 
creases. 

The  data  of  the  present  paper  enable  no  definite  conclusions 
to  be  drawn  regarding  the  validity  of  the  drag  correction. 
However,  in  view  of  the  accuracy  of  the  other  corrections  for 
the  circular  tunnel  and  in  view  of  the  fact  that  tlu'  corr(>- 
sponding  drag  correction  for  a  two  dimensional  tuumd  is 
known  to  be  accurate,  it  is  to  be  expected  that  this  correction 
will  give  a  satisfactory  evaluation  of  the  wall  interference 
upon  the  measured  drag. 

The  equations  of  the  present  jiaper  should  not  be  expected 
to  give  accurate  residts  wlien  applied  to  tests  in  which  air 
leakage  occurs  at  the  tunnel  walls.  In  such  tests  the  lift  at 
the  walls  drops  markedly,  so  that  the  assumption  that  the 
lift  is  uniform  across  the  span  is  no  longer  valid.  The 
importance  of  avoiding  such  leakage,  if  reliable  airfoil 
characteristics  are  to  be  obtained,  is  pointed  out  in  reference 


5  with  regard  to  tests  in  two-dimensional  tunnels.  The 
same  general  considerations  apply  in  the  case  of  an  airfoil 
spanning  a  circular  tunnel. 

CONCLUSIONS 

Airfoil  data  obtained  from  tests  at  subsonic  speeds  of  an 
airfoil  spanning  the  center  of  a  closed-throat  circular  wind 
tunnel  can  be  coirected  to  free-air  conditions  by  means  of 
the  following  equations; 


t— I  [1  —  1— (Af')^ 

tl- </  (  '  Al-(Af')- 


(48) 


/.  O.  (  ,  ,  I  ->A7  (Ary  [1  -  (),7(A/')d[l  1-0.4(A/')^]  ,) 


(49) 
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where  tj,  <ti,  and  a-  are  given  by 

(j) 


(1.289  (51 


cr^  —  0 


-(0 


(46) 
(55) 

(47) 


and  A  is  a  dimensionless  factor  the  value  of  which  depends 
upon  the  shape  of  the  base  profile  of  the  airfoil.  (See  table 
I  and  equation  (3)  of  reference  5.)  The  remaining  symbols 
arc  defined  in  Appendix  C.  Numerical  values  of  the  func¬ 
tions  of  M'  which  appear  in  these  equations  are  given  in 
table  II.  Experimental  pressure  distributions  can  also  be 
corrected  by  proper  modification  of  the  method  of  reference  5 
as  indicated  in  the  text. 

Tests  of  an  N.AC.A  4412  airfoil  at  low  speed  for  two  ratios 
of  airfoil  chord  to  tunnel  diameter  demonstrate  the  valiflity 
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of  the  foregoiog  eqiiatioDS  at  low  Mach  numbers.  At  angles 
of  attack  below  the  region  of  maximum  lift,  the  equations  are 
applicable  for  chord-diameter  ratios  up  to  at  least  0.625,  the 
maximum  ratio  tested.  In  the  region  of  maximum  lift  a 
chord-diameter  ratio  of  0.35  is  known  to  be  permissible,  and 
still  higher  ratios  may  give  satisfactory  results.  An  examin¬ 
ation  of  the  validity  of  the  equations  at  high  Mach  numbers 
is  not  possible  at  present,  but  the  maximum  permissible 
chord-diameter  ratios  may  be  expected  to  decrease  as  the 
Mach  number  increases. 

The  tests  also  indicate  that  at  low  Mach  numbers  the  span- 
wise  lift  distribution  on  an  airfoil  spanning  a  closed-throat 
circidar  tunnel  is  essentially  constant  except  at  angles  of 
attack  in  the  immediate  vicinity  of  the  stall.  This  result 
corroborates  the  general  conclusion  of  Appendix  B,  in  which 


it  is  demonstrated  that  the  lift  is  uniform  across  an  untwisted, 
constant-chord  airfoil  spanning  any  closed-throat  wind  tun¬ 
nel,  irrespective  of  the  cross-sectional  shape  of  the  tunnel. 

The  correction  equations  cannot  bo  expeiited  to  apply  at 
or  in  the  immediate  vicinity  of  the  choking  Mach  number, 
which  is  the  maximum  Mach  number  attaii'able  with  a  given 
combination  of  airfoil  and  tunnel  test  section.  The  choking 
Mach  number  can  be  estimated  by  means  of  equations  given 
in  the  report. 
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APPENDIX  A 


TRANSFORMATION  OF  SERIES  OF  BESSEL  FUNCTIONS 

Th  0  soi'ies  involving  Bessel  functions  which  appear  in  the 
discussions  of  the  interference  effects  associated  with  airfoil 
•  camber  and  thickness  are,  .as  pointed  out  in  the  text,  pooi'ly 
suited  for  use  at  small  values  of  the  variable  x.  It  will  be 
shown  here,  by  means  of  a  method  demonstrated  by  Watson 
(reference  15),  that  the  series  maj^  each  be  expressed  as  a 
combination  of  elementary  functions  and  a  convergent  power 
series.  The  resulting  scries  arc  well  adapted  for  use  in  the 
present  problem.  The  notation  used  for  the  Bessel  functions 
is  that  of  Watson  (reference  Id)  and  of  the  Smithsonian 
Tables  (reference  14). 

Series  for  camber  effect. — The  discussion  of  the  interfer¬ 
ence  effects  associated  with  airfoil  camber  involves  the  series 


HV 


.^^X,,/•Jl(X,,r),//'(Xsr) 


(AD 


corn'crgent  for  negative'  values  of  .r.  The  sumnialioii  with 
respect  to  k  extends  over  all  the  positive  roots  of  the  (‘({uatioji 

J/(X,,r)=0  (A2) 


Betting  ,/5  =  X,,/'  and  k~~x,  tlu'  se'ries  may  be  w  iitleii 
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(Ad) 


wliere  the  summation  is  taken  over  all  tlie  positive  roots  of  tlie 
c‘r|uatiou 

</i'(.h)-0  (A4) 


Now,  consider  the'  function 

_£  J\'  iw)Y\  (n-r]lr)  —  Ji{'ii-r]jr)  B/BD 

2  Ji\w) 


{Ary) 


where  the  quantity  I’l  is  a  Be'ssel  function  of  tlu'  second  kind 
of  order  miit}'.  This  function  has  a  simple  pole  at  <‘ach  of 
the  ])oints  w=  ±js  and  is  one-valued  and  analytic  at  all  other 
points  in  tlu'  com])le.\  w-])lane.  Its  :residue  at  the  point  j, 
can  be  shown  to  be 


Jii.jsvlr)  ^ 

isJ 1  (js)Jl"  ijs) 


which  is  identical  with  the  general  tei’in  of  tlie  series  (A3). 
By  th('  theorem  of  residues,  the  integral  of  the  function  (A5) 
taken  counterclockwise  around  a  contour  inclosing  the  por¬ 
tion  of  the  complex  plane  to  the  right  of  the  imaginary  axis 
is  then  equal  to  27rflT).  The  integral  along  a  large  semi¬ 
circle  on  the  right  of  the  imaginary  axis  tc'iids  to  zero  when 
the  radius  of  the  semicircle,  tends  to  infinity  through  values 
such  that  the  semicircle  avoids  the  poles  of  the  integrand. 
It  is  thus  necessary  to  retain  only  the  integral  along  the 
imaginary  axis.  Tlu'  contour  must,  however,  have  an  in¬ 


dentation  to  the  right  of  the  origin,  since  the  integrand  has 
a  pole  there  with  residin'  (r^  +  7j^)/rj7.  If  the  radius  of  the 
indentation  is  made  to  approach  zero,  ITj  may  finally  be 
written 


1^1  = 


2rr] 


Jj{w)  Yy  (wnlr)—Jiiwijlr)Yi'{w) 


e-^’^^rd-w 


(AO) 


From  the  known  relations  for  the  modified  Bessel  functions, 
it  is  jcadily  shown  that 


1 1(±'(0  =  — '-fi(0±'"  ^A|(0 

TT 

<^d(±  /()  =  //(() 

1  /(± 'd)=  ± ?// (0  A”  AV (() 

TT 


(A7) 


wliei-e  /i  and  Ky  are  modified  Bessel  functions  of  the  first 
and  sc'cond  kind  of  order  unity.  By  writing  the  integral 
in  equation  (A6)  in  two  parts,  one  along  the  positive  and 
OIK'  along  tlu'  iK'gaiive  imaginary  axis,  and  rc'placing  w;  in 
tiu'se  inh'grals  by  +it  and  —if,  respectively,  Wy  then 
liecomes 


ir, 

or 


I  1  f'”  Y'ii)Ky(Jr}lr)  —  I,{tri/r)Ky'{t) 
2rri  ir  J,,  //(() 


sin  {tKjr)dt 


U,- 


Ky{trilr)  sin  (t/c/r)  di— 


Ky'it) 

Y'{t) 


ly  (tr]/r)  sin  (N'/r)  dt 


(AS) 


The  valiu'  of  the  first  inti'gral  in  this  ('quation  is  given  by 
Watson  in  reference  15  as 

[  Kyitvlr)  sin  (tn/r)  dt  =  ^  ,,  (A9) 


Tli(^  second  integral  can  be  evaluated  by  expanding  the 
proiluct  lyitrijr)  sin  (t/c/r)  in  ascending  powers  of  t  and  inte¬ 
grating  term  by  term.  The  series  expansion  for  the  prod¬ 
uct  is 

Iy{iv/r)  sm  {tK/r)  ^  j.,  (2p+Yl  2-" '  ‘  '  i) 

and  the  term-by-term  integration  gives 


Ky\t) 

h'W 
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Yitv/f)  sin  {tK/r)  dt 
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The  coeflicient  =  's  given  l)y 


,  f”  /v,'(/) 

jo 


M  2/ 


(0 


t-'  (li 


wliieli  may  be  written  after  integration  by  ]iarts 

(2/+l)2rj,, 

This  integral  is  a  constant  for  any  given  value  of  f. 

Reverting  to  the  oi'iginal  variable  ,r,  tbe  e.xpansion  for  ll'i 
may  finally  be  written 

‘  2r?)  277  + 


f.  (“i)'’  /2'2(*+„+,)  r‘  " 
l-\  (*+  1 ) !  (2p+  1 ) !  2^‘-+' 


(A  12) 


This  agrees  witli  tbe  result  given  without  derivation  by 
Tani  and  Taima  (reference  18). 

For  purposes  of  computation  the  coefficient  /x'2/  is  writ  ten 


vberc 
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(Aid) 


(A  14) 


1'bc  ciuantity  &'2f  can  then  be  e.xprossed  in  a  form  suitable 
for  computation  by  means  of  a  method  devised  by  Watson 
for  an  analagous  integral  (reference  15). 

As  the  first  step,  tbe  function 

ff 

[//(OPcosfTTt/f) 


is  written  as  a  sum  of  partial  fi'actions,  h  being  it  positive 
constant  which  will  be  fixed  later.  Fbis  can  be  accomplished 
by  considering  the  integral 


r  _  vfUhv 

j  (w— 0  {Ii{w)\ -  cos  (TrM-’/6) 


(A  10) 


tiround  the  circle  \w\—li  in  the  complex  ]dane.  The  inte- 
gi-and  of  the  integral  (AlO)  has  jioles  at  the  points 

W=f,  w=±  ijs,  W=  ±  («+ 1  /2)b 


where  j,  is  a  positive  zero  of  J^'iw);  and  8=1,  2,  d,  ..., 
«  =  0,  1,  2,  d,  .  .  .  .  The  residue  at  the  simph^  jiole  at  w-  ( 

is  the  function  (j\.15).  The  residues  at  the  sim[)l('  jroles  iit 
1/’=  ±(?(,+  l/2)6  are 

(- I )''■(«■+ 1/2)''  ()'•'+' 

7r(n6+  6/2  =Ft)  [l7\nb  +  h/2  jY 


The  i)oh's  at  i/;=  i vjs  ai’c  second  oi'dcr  poles;  tlu^  residues 
there  are 

(—1  r 

J7{js)  (1  — j/)'  cosh  (7rj,/h)  (./»±  d)' 

Zf-i-^js/b)  tanh  (7rj,,/6)  +  '^  j". 


Now,  the  integral  (AlO)  taken  around  the  circle  lw|  =  /f 
tends  to  zero  when  E  tends  to  iuliuity  in  such  a  maimer 
that  the  circle  never  passes  through  a  pole  of  the  iutogrand. 
It  follows  from  the  theorem  of  residues  that  the  sum  of  the 
residues  of  the  integrand  at  all  its  poles  is  zero;  thus 


ff  2  ^  _  (_^1)^  (ft+l/2)''+'  6'-'+' 

cos^(7rt/6)"=  X  h  [(/i6  +  6/2)'-F]  [//(a6+6/2)l' 
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fly  multiplying  tins  ('quation  by  cos  (Tr</6)  and  inti'grating 
fiom  —  oD  to  +  00^  it  can  be  shown  with  the  aid  of  certain 
integral  relations  given  by  Watson  (reference  15,  p.  dfj)  that 


m  ^  (n  +  ]/2)'^6'w  . 

J_.  [//(/)F  ^^,[7/(n6+6/2)F 

CO  j  — 
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and  therefore 


1  f”  __6  ^  {,nh  +  bl2)-f 

^  [imY~^h  W{nh  +  hl2)Y^ 

JFliJ(T-J?)'  [coiS(S/&r 
(2/-+  \  J-^)  J  ' '' ^ 

'riie  first  sei-ies  in  this  equation  converges  rapidly  when  b  is 
large,  the  second  when  b  is  small.  A  reasonable  compromise 
for  purposes  of  calcuation  is  to  take  6=1. 

Equation  (AlO)  with  6=1  has  been  used  together  with 
equation  (Aid)  to  determine  the  first  four  vahu's  of  the 
coefficient  m'>/.  The  final  results  arc 


-0.999 
Ij.\=---1A)27 
72^6=  —  9.78 
72'8=- 120.8 


(A2()) 


Gomparable  values  of  n'i  and  ij.\  to  the  same  number  of 
significant  figures  are  given  witfiout  derivation  in  reference 
18.  The  value  of  /u'2in  this  latter  reference  agrees  wdth  that 
of  the  present  paper  but  /rh  differs  by  one  in  the  third  di'cimal 
j)lacc.  The  value  given  in  (A20)  has  been  carefully  checked 
for  several  values  of  the  parameter  6  and  appears  to  be  cor¬ 
rect.  Values  of  g'e  and  m's  apparently  bave  not  previously 
been  computed. 

Series  for  thickness  effect. — The  series  which  appears  in 
the  discussion  of  the  interference  effects  associated  with 
airfoil  thickness  is 


_  '^F*5/,(X.,77) 

d„'(X,/-) 


(A21) 
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convergent  for  negative  values  of  x.  The  summation  with 
respect  to  s  extends  overall  the  positive  roots  of  the  equation 


Reverting  to  the  original  variable  x,  the  expansion  for  H'o 
may  finally  be  written 


Ji(X,,r)  =  0  (A22; 

Lett  ing  j,<=X,r  and  k=  —x  as  before,  the  series  may  he  writt  en 


II  "  \ 

^  Zj  -/oU)  ' 


-i.n/r 


iA2'A) 


where  the  summation  is  takmi  oviU'  the  ))osilive  Tools  of 

-  0  (A24) 


'rile  function 

TT  J ,  (w)  Y 1  (U’rj/r)  —  J ^ {wn/r)  Ti ( w) 

2  J,{w) 


(A2r,, 


lias  a  simple  pole  at  each  of  the  points  w=  ±js.  Its  residue 
at  each  of  these  points  can  he  shown  to  he  identical  with  the 
gimeral  term  of  the  series  (A2d).  Unlike  the  function  in  the 
previous  series,  this  function  is  regular  at  the  origin.  Inte¬ 
gration  around  the  portion  of  the  complex  plane  to  thi'  right 
of  till'  imaginary  axis  then  gives 


ir, 


2x1  J  _ 


'  X  Jiiw)  Yi  (ii'ii/c)  I'l  (w) 

xi  2  JAw) 


i/’e  “  "''••Ulw 

(A2r.) 


By  applying  the  first  two  of  (>quations  (A?)  and  combining 
the  inti'grals  along  tlu'  two  lialvi's  of  the  imaginary  axis  as 
liefori',  the  series  becomes 

1  1'“  1  I’"  K  H'\ 

ll'i  -  tKAtvIA  cos  {hlr)(lt~  -  fliUnli')  cos  (lujAdt 

irjii  xjd  iitt) 

(A27) 


The  first  integral  can  he  evaluatial  by  difl'erentiating  rela¬ 
tion  (A9)  witli  respect  to  k.  'Pliis  oiicration  gives 

COS  =  (A2S) 


'I'he  second  integral  can  be  evaluated  as  before  by  expand¬ 
ing  the  product  tliAvh')  cos(L//')  in  ascending  iiowers  of  t 
and  integrating  tei-m  by  term.  'Lhe  series  expansion  for  the 
product  is 


CO  CO  ( _ 1  Ip /2(*+r+n„2A-+i 

tlSnA')  cos  (A/c)  ]),  (-2^)! 2-'*  + 

and  the  term-by-term  integration  gives 

li  {tK/r)di 


where  t  he  coefficient  p4u  =  M2/  is  given  by 

1 


^2/ 


,  p 


K\{i)  ,;i  I, 

7,(0  ■(2/+I) 


lAdf 

T^J»  lA 


it) 


(AdO) 


The  integral  (AdO)  has  been  investigated  by  Watson  (re¬ 
ference  15).  Its  value  for  any  given  r  can  be  computed  from 
the  series 


(271  Dmv 


d'dt 
At) 


1  p  /2, 

Id 

_CC^  12/ -I  p 

^  'S-LTjd  [e 

{21+ A 


Aih  +  bi'!)-'' 
irf+u  !,-{>+  \  h/-2]~' 


+jAt>) 

osh"(xy,//>) 


1 


cosh  +3  Ah)} 


where  h  is  an  arbitrary  iiositlve  constant.  'Phis  equation  has 
heim  used  with  6  -  1  to  didermine  the  first  four  values  of 
'Pile  final  results  are 


0.797 

Ml -=1.200 

M6=7.4fi 
Ms  =90.2 


(Add) 


'Pile  (irst  two  of  these  vahu's  agri'c  to  the  three' decimal  ])lae('s 
with  the-  two  numerical  values  computed  by  Watson.  'Pin' 
ri'maining  two  vahu's  lia-vi'  not  jn'eviously  been  eomjmted. 


APPENDIX  B 

CONSTANCY  OK  LIFT  OVER  AN  AIRFOIL  SPANNING  A  CLOSED- 
THROAT  TUNNEL 

Consider  an  infinitesimally  thin  untwisted  airfoil  of  con¬ 
stant  choi'il  spanning  a  closed-throat  wind  tunnel  of  arbitrary 
section.  Such  an  arrangemi'iit  is  shown  in  figure  10,  which 
is  a  section  of  the  tunnel  as  seen  from  downstream.  It  is 
assumed  that  the  flow  in  the  tunnel  is  nonviscous  and  tbat 
the  airfoil  therc'foie  has  no  drag. 

Suppose  for  the  time  being  that  the  lift  varic's  in  some 
manner  across  the  span  of  the  airfoil.  Any  such  valuation 
will  be  accompanied  by  a  system  of  vortices  trailing  from  the 
ail-foil  and  extending  infinitely  far  downstream.  If  the  usual 
assumption  is  made  that  the  trailing  x"ortiees  ari'  parallel  to 
the  axis  of  the  tunnel,  the  flow  pattei-n  in  a  plane  normal  to 
the  axis  at  infinity  downstream  must  lie  of  the  nature  shoivn 
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in  figui’e  10.  The  flow  pattern,  in  general,  consists  of  a  nuTii- 
ber  of  separate  sections  within  each  of  which  the  flow  has  a 
closed,  circulatory  charactei'.  The  line  AB,  which  repj'esents 
the  projection  of  the  airfoil,  extends  across  every  such  section, 
and  each  of  the  sections  contains  the  filaments  of  a  portion  of 
the  system  of  trailing  vortices.  The  exact  character  of  the 
flow  pattern  in  any  particular  case  depends  upon  the  span- 
wise  variation  in  lift  and  upon  the  cross-sectional  shape  of  (he 
tunnel. 

Now,  consider  the  flow  around  a  streamline  within  any  one 
of  the  separate  sections  of  the  flow  pattern — say  the  stream¬ 
line  CD  in  the  section  at  the  left-hand  side  of  the  tunnel  in 
ligure,  10.  This  streamline,  like  all  the  sti-eamlines,  inter- 
sc'cts  the  projection  AB  of  the  airfoil  in  two  points,  denotrul 
as  C  and  D  in  the  figure.  The  fact  that  in  the  presence  of 
the  tunnel  walls  each  streamline  must  intersect  AB  in  two 
])oints  is  essential  to  the  discussion.  If  it  is  supposed  for 
purposes  of  discussion  that  the  direction  of  flow  is  clockwise 
as  indicated,  the  vertical  component  of  velocity  at  C  is 
upward  while  the  corresponding  component  at  D  is  down¬ 
ward.  This  direction  of  flow  corrc'sponds  to  a  net  circulation 
in  the  clockwise  direction  for  all  the  trailing  vortex  filanuuits 
enclosed  within  the  streamline. 

At  the  position  of  the  airfoil  the  pattern  of  transvei-se 
velocities  induced  by  the  trailing  vortices  is  gt'ometiically 
similar  to  the  pattern  at  infinity  dowtistream,  only  tin? 
magnitude  of  the  velocities  being  different.  Hence,  at  points 
on  the  airfoil  directly  ahead  of  point  C,  the  vertical  velocity 
induced  by  the  ti'ailing  vortices  is  upward.  At  ])oints 
directly  ahead  of  point  D,  the  velocity  is  downward.  Tims, 
since  the  airfoil  is  untwisted,  the  airfoil  section  corresponding 
to  C  operates  at  a  largc'r  efi'ective  angh;  of  attack  than  does 
the  section  corresponding  to  D.  If  the  aiifoil  is  of  constant 
chord  as  assumed,  this  means  that  the  lift  at  section  C  jnust 
he  greater  than  the  lift  at  section  D. 

As  has  been  pointed  out,  however,  the  trailing  vortices 
discharged  between  sections  C  and  D  must  lia\ c  a  net  cii’cu- 
lation  in  the  clockwise  direction  iti  figure  10.  This  means 
that  the  circulation  of  the  spanwise  hound  vortices  at  sec¬ 
tion  D  must  be  greater  than  at  section  C.  Since  the  direc¬ 
tion  of  stream  flow  was  taken  to  be  toward  tlu'  obs('i'ver,  this 
in  turn  means  that  the  lift  at  section  C  must  be  less  than 
that  at  section  D,  which  is  in  direct  contradiction  to  the  pre¬ 
vious  resrdt.  The  original  supposition  that  the  lift  varies 
across  the  span  thus  leads  to  two  mutually  contradictory 
conclusions  and  is  therefore  invalid.  It  follows  that  the 
spanwise  distribution  of  lift  is  uniform  across  an  untwisted, 
constant-chord  airfoil  spanning  any  closed-throat  wind 
tunnel,  irrespective  of  the  cross-sectional  shape  of  tlie 
tunnel. 

As  mentioned  at  the  outset,  this  result  depends  upon  tlu' 
assumption  that  the  airfoil  is  infinitesimally  thin  and  has  no 
drag.  It  will  not  be  strictly  true  if  tlu'  inci’case  in  ('ffective 
stream  velocity  caused  by  the  interference  between  the  walls 
and  the  airfoil  thickness  and  wake  is  not  uniform  across 
the  span.  The  result  also  neglects  any  effect  that  the  bound¬ 
ary  layer  along  the  walls  of  the  tunnel  may  have  upon  the 
lift  distribution.  That  these  approximations  are  not  serious, 
at  least  in  the  case  of  the  circular  tunnel,  is  indicated  by  the 
experimental  results  of  figures  6  and  7. 


The  foregoing  reasoning  is,  of  course,  inapplicable  for  an 
airfoil  which  does  not  span  the  tunnel  or  for  a  finite-span 
airfoil  in  free  ah-.  In  thesi'  instances,  the  pi'ojection  of  the 
airfoil  docs  not  extend  across  all  of  the  sections  into  which 
tile  transverse  flow  pattern  is  divided,  and  the  streamlines 
of  this  pattern  need  not  intersect  the  projection  of  the  airfoil 
in  twm  points.  Under  tlu'se  conditions  a  type  of  varying  lift 
distribution  can  he  found  which  does  not  lead  to  a  logical 
inconsistency. 

APPENDIX  C 

LIST  OF  IMPORTANT  SYMBOLS 

c  airfoil  chord 

</  diameter  of  circular  tunnel 

/■  i-atlius  of  circular  tunnel 

li  height  of  rectangular  tunnel 

/u  height  of  rectangular  tunnel  equivalent  to  circular 

tunnel  with  regard  to  camber  effect 
li,  height  of  rectangular  tunnel  equivalent  to  circular 
tunnel  with  regard  to  thickness  effect 

T  1/4  chord-height  factor  wdth  regard  to  wake 

eft'ect  in  rectangular  tunnel 

T>  0.421  chord-diameter  factor  with  regard  towndee 
effect  in  circular  tunnel 

<T  ^  ^  Aliord-height  factor  with  regard  to  camber 

and  thickness  effect  in  rectangular  tunnel 
<T|  0.289  ;  chord-diameter  factor  with  regard  to 

camber  efl’ect  in  circular  tunnel 
IT-,  0.449  ;  chord-diameter  factor  with  regard  to 

thickness  effect  in  circular  tunnel 
factor  depending  upon  shape  of  base  profile  (see 
equation  (44)  and  table  I) 
a  angle  of  attack 

Ci  section  lift  coefficient 

section  quarter-chord-moment  coefficient 
c,,  section  drag  (oefficient 

U  stream  velocity 

<1  dynamic  pressure 

l\f  Mach  number 

li  Reynolds  number 

.r,y,z  I'ectangular  sjiace  coordinates 

cylindrical  space  coordinates  (see  equations  (1)) 
r'  circulation  of  single  line  vortex  in  tunnel 
ilV'ld^  circulation  per  unit  chord  length 
)7,f  y  and  z  coordinates  of  elementary  vortex 
coii,0,i  radial  and  angular  coordinates  of  elementary  voidex 

$  chordwisc  coordinate  of  elementary  vortex;  also 

variable  of  integration  in  equations  (2)  and  (4) 

0  velocity  potential 

r/ ,r/  X  and  s  components  of  induced  velocity 
Ah''  increase  in  axial  velocity  at  position  of  airfoil  in 
tunnel 

IX  doublet  strength 

1-1,  projected  thickness  of  airfoil 

effective  thickness  of  airfoil 
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irijH'i)  series  ol  terms  involviiifr  Bessel  functions  (see  equa¬ 
tions  (At)  and  (A21)  of  Appendix  A) 

Jm,yn  Bessel  functions  of  first  and  second  kind  of  oi'der  m 
(Watsonk  notation) 

Iiii.Km  modified  Bessel  functions  of  first  and  second  kind  of 
order  m  (Watson’s  notation) 

X,5  variable  of  summation  defined  by  tlie  roots  of  the 
e(|uation  t/,/(Xs/')  =  0 

js  X,,/-;  root  of  the  einiation  J,,/ (,7,)  =  d 

„j.  numerical  coefficients  (see  equations  (doj  and  (13)) 
variables  of  summation 
t,w  variables  of  integ-ration 
K  alternate  variable  defined  as  e(|ual  to  .r 

Siqierscripts 

(')  wIhui  jiertaining  to  fluid  pro])erties,  denotes  values 
in  the  undisturhed  stream  in  the  tunnel;  when 
pertaining  to  airfoil  characteristics,  denotes  values 
in  tunnel,  coefficients  being  refei'red  l.o  dynamic 
])ressure  q' \  denotes  fii'st  derivative  of  Bessel 
function  with  respect  to  its  argument 
(")  denotes  second  derivative  of  Bessel  function  with 
respect  to  its  argumeid, 
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Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

LonsfitiidiMal 

X 

X 

Rolling _ 

L 

Y - *Z 

Roll 

u 

P 

Lateral . 

Y 

Y 

Pitching 

M 

Z - *X 

Pitch . 

e 

V 

? 

Normal 

Z 

z 

Vftwing 

N 

X - ►F 

Yaw 

'P 

ID 

r 

Absolute  coefficients  of  moment 

p  _ p _ 

'’‘~qcS 

(rolling)  (pitching)  (yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5,  (Indicate  surface  by  proper  subscript.) 


4,  PROPELLER  SYMBOLS 


D 

P 

Diameter 

Geometric  pitch 

P 

p 

Power,  absolute  coefficient  pj^zjyi 

pID 

v 

Pitch  ratio 

Inflow  velocity 

c. 

Speed-power  coefficient  = 

Slipstream  velocity 

V 

Efficiency 

T 

Q 

—  •  T 

Thrust,  absolute  coefficient  Cr— — 

pn‘D^ 

Torque,  absolute  coefficient 

n 

4- 

Revolutions  per  second,  rps 

Effective  helix  anglo=tair‘^2“^^^ 

1  hp=76.04  kg-m/s=5.50  ft-lb/aec 
1  metric  horsepower=  0.9863  hp 
1  mph=0.4470  mps 
1  mps= 2.2369  mph 


5.  NUMERICAL  RELATIONS 

1  lb=:0.4536  kg 
1  kg=2.2046  lb 
1  mi=  1,609.35  m= 5,280  ft 
1  m= 3.2808  ft 


The  Effect  of  Wall  Interference  upon  the  Aerodynamic  Characteristtcs  of  an  Air¬ 
foil  Spanning  a  Closed  Throat  Circular  Wind  Tunnel 
>R(S):  Vlncentl.  Walter  G.;  Graham,  Donald  3. 
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